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Summary. These lecture notes address a central problem of theoretical nuclear 
physics; how to establish a relationship between low-energy, non-perturbative QCD 
and nuclear phenomenology which includes both nuclear matter and finite nuclei. 
We develop a microscopic covariant description of nuclear many-body dynamics 
constrained by chiral symmetry and in-medium QCD sum rules. A relativistic point- 
coupling model is derived, based on an effective Lagrangian with density-dependent 
contact interactions between nucleons. These interactions are constructed from chi- 
ral one- and two-pion exchange, combined with the large isoscalar nucleon self- 
energies that arise through changes in the quark condensate and the quark density 
at finite baryon density. Nuclear binding and saturation are almost completely gen- 
erated by chiral (two-pion exchange) fluctuations in combination with Pauli effects, 
whereas strong scalar and vector fields of about equal magnitude and opposite sign, 
induced by changes of the QCD vacuum in the presence of baryonic matter, gen- 
erate the large effective spin-orbit potential in finite nuclei. Promising results are 
found for the nuclear matter equation of state and for the bulk and single-nucleon 
properties of finite nuclei. 



1 Introduction 

The physics based on Quantum Chromodynamics (QCD), the fundamental theory 
of interacting quarks and gluons, is undoubtedly one of the most fascinating fields 
of modern science. It covers an enormously wide spectrum of phenomena, ranging 
from hadrons and nulcei to matter under extreme conditions of temperature and 
density. How the basic quark constituents and gluon fields arrange themselves to 
form nucleons and mesons, and how these localized clusters of confined quarks and 
gluons interact collectively in nuclei: these are key questions of nuclear physics. 

At very short distance scales, r < 0.1 fm, corresponding to momentum transfers 
above several GeV/c which probe space-time intervals deep inside the nucleon, 
QCD is a theory of weakly interacting, pointlike quarks and gluons. The rules 
for their dynamics are set by local gauge invariance under SU{3)coior- This is the 
area of perturbative QCD, an active field of research all by itself but with little 
direct impact on the understanding of nuclei in which the relevant length scale 
(e.g. the average separation between nucleons) is typically more than twice the size 
of the individual nucleon. At such large distance scales, r > 1 fm, QCD is realised 
as a theory of pions coupled to nucleons (and possibly other heavy, almost static 



2 Dario Vretenar and Wolfram Weise 

hadrons). Their dynamics is governed by an approximate symmetry of QCD: Chiral 
Symmetry, and its spontaneous breakdown at low energy. 

The non-linear dynamics of gluons and their strong interactions with quarks 
imply a complex structure of the QCD ground state (or "vacuum"): it hosts strong 
condensates of quark-antiquark pairs and gluons, an entirely non-perturbative phe- 
nomenon. The quark condensate (qq), i.e. the ground state expectation value of 
the scalar quark density, plays a particularly important role as an order parameter 
of spontaneously broken chiral symmetry. Hadrons, as well as nuclei, are excita- 
tions built on this condensed ground state. Changes of the condensate structure 
of the QCD vacuum in the presence of baryonic matter are a source of strong 
fields experienced by the nucleons. The nuclear dynamics governed by these strong 
fields, together with the pion-exchange forces that determine the nucleon-nucleon 
interaction at long and intermediate distances, will be important elements in our 
discussion. 

Our framework is the effective field theory that represents low-energy QCD with 
the lightest {u and d) quarks. This theory is written in terms of an effective La- 
grangian constructed according to the symmetries and symmetry breaking patterns 
of QCD. Identifying the relevant, active degrees of freedom at low energy is the first 
important step in this approach. Nuclei are aggregates living in the low-temperature 
(hadronic) sector of the QCD phase diagram, where quarks and gluons are confined 
in color singlet mesons and baryons. At the energy and momentum scales charac- 
teristic of nuclei, only the states of lowest mass in the meson and baryon spectrum 
are expected to be "active": pions and nucleons. Their quark-gluon substructure 
and detailed short-distance dynamics are not resolved at nuclear scales. In the ef- 
fective field theory the substructures of pions and nucleons are encoded through 
measured structure constants and through contact interactions with parameters to 
be determined empirically. The effective field theory approach does have predictive 
power, in the sense that there exists a guiding (symmetry) principle for construct- 
ing an ordered hierarchy of terms in the effective Lagrangian. The parameters of 
the leading terms are controlled or constrained by measurements of independent 
subprocesses such as pion-nucleon scattering. 

Nuclear phenomenology has, of course, a long history and a rich empirical back- 
ground that we can build on. Concepts such as density functional theory and effec- 
tive field theory are at the basis of many successful nuclear structure models. Quan- 
tum Hadrodynamics (QHD), in particular, represents a class of Lorentz-covariant, 
meson-nucleon or point-coupling models of nuclear dynamics, that have been used 
successfully in studies of a variety of nuclear phenomena 1. . Models based on the 
relativistic mean-field (RMF) approximation of QHD have been applied to describe 
various properties of spherical and deformed nuclei all over the periodic table 

Pions, as Goldstone bosons of spontaneously broken chiral symmetry, are well 
known to be important ingredients in nuclear systems The behavior of the 
nucleon-nucleon interaction at long and intermediate distances is determined by 
one- and two-pion exchange processes. Ab-initio many-body calculations of nu- 
clear matter and light nuclei ^ strongly emphasize the crucial role played by the 
pion-exchange tensor force. Nevertheless, pions are usually not included as explicit 
degrees of freedom in QHD-type calculations. The most successful QHD models 
are so far purely phenomenological. Effects of two-pion exchange are treated only 
implicitly through an ad-hoc scalar-isoscalar mean field. 
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The successes of relativistic nuclear mean-field approaches have been attributed 
primarily to the large Lorentz scalar and four- vector nucleon self-energies which are 
at the basis of all QHD models [S]- There is substantial (though indirect) evidence, 
in particular from nuclear matter saturation and abnormally large spin-orbit split- 
tings in nuclei, that the magnitude of each of the scalar and vector mean fields 
experienced by nucleons in the nuclear interior is of the order of several hundred 
MeV. Studies based on finite-density QCD sum rules have shown how such large 
scalar and vector nucleon self-energies can arise naturally from the in-medium 
changes of the scalar quark condensate and the quark density |^ |7| In most 
applications of QHD models to finite nuclei, however, the interactions generating 
these large nucleon self-energies have been treated in a purely phenomenological 
way. In some models they arise from the exchange of "sigma" and "omega" bosons 
with adjustable coupling strengths. The physical meaning of such assumed ad-hoc 
mechanisms may be seen in different perspective when exploring their basic QCD 
connection. 

These lecture notes describe and summarize recent developments based on the 
following conjectures that are believed to establish links between low-energy QCD, 
its symmetry breaking pattern, and the nuclear many-body problem: 

- The nuclear ground state is characterized by large scalar and vector fields of 
approximately equal magnitude and opposite sign. These fields have their origin 
in the in-medium changes of the scalar quark condensate (the chiral condensate) 
and of the quark density. 

- Nuclear binding and saturation arise predominantly from chiral (pionic) fiuctu- 
ations (reminiscent of van der Waals forces) in combination with Pauli blocking 
effects, superimposed on the condensate background fields and calculated ac- 
cording to the rules of in-medium chiral perturbation theory. 

We begin by briefiy summarizing important aspects of low-energy QCD and chiral 
symmetry inasmuch as they are relevant to our theme. Elements of in-medium 
chiral perturbation theory and results from QCD sum rules at finite density are also 
described. We then proceed by studying nuclear matter and nucleon self-energies in 
the context of chiral dynamics and with constraints from the expected changes of 
condensates at finite baryon density. In the next step, these features are translated 
into an equivalent relativistic point-coupling model with density dependent coupling 
strengths. This model is used to calculate a variety of bulk and single particle 
properties of finite nuclei, from ^^O to ^''Zr. We conclude with an outlook and a 
discussion of further steps that are still ahead. 

Much of the material and many of the results presented here have appeared in 
recent publications |^ llUl[mil2|[T^ . A presentation of the general framework can 
be found in previous lecture notes 1141 . Introductory surveys on basic QCD notions 
related to our subject are given in several chapters of ref. 1151 . 
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2 Low-energy QCD 
2.1 Chiral symmetry 

Our framework is QCD in the sector of the lightest {u, d) quarks. They form a 
flavour Nf = 2 (isospin) doublet with "bare" quark masses of less than 10 MeV. 
The flavour (and colour) components of the quarks are collected in the Dirac fields 
il>{x) = {u{x) , d{x))'^ . The QCD Lagrangian, 

Cqcd = ^(iTM^" - - ^TriG^^.G'^n, (1) 

involves the SU{3)coior gauge covariant derivative D'^ and the gluonic field tensor 
= {i/g)[D'^ , D"]. The 2x2 matrix m — diag{mu, rud) contains the light quark 
masses. The strange quark is more than an order of magnitude heavier (rris ~ 150 
MeV), but still sometimes considered "light" on the typical GeV scales of strong 
interaction physics. It will nevertheless be ignored in the present discussion, together 
with the heavy quarks {Q — c, b and t). 

Consider QCD in the limit of massless quarks, setting m = in Eq.lQ. In this 
limit, the QCD Lagrangian has a global symmetry related to the conserved right- 
or left-handedness (chirality) of zero mass spin 1 /2 particles. Introducing right- and 
left-handed quark fields, 

V'fl.t = ^(1±75)V', (2) 
we observe that separate global unitary transformations 

V'fl exp[z6l^y] tpR, tpL exp[i6l£y] ^pL, (3) 

with Ta{a — 1,2,3) the generators of (isospin) SU{2), leave Cqcd invariant in the 
limit m ^ 0. This is the chiral SU{2)r x SU{2)l symmetry of QCD. It implies 
six conserved Noether currents, Jr ^ ~ ^Rl^^-'^R and a, ~ ''Pl'y'^ -^ipL, with 
dfiJ^ = dfj^J'^ = 0. It is common to introduce the vector and axial currents, 

Va^ = J'ia + Jla = V'l" f , ^C^^) = J R,a - JL,a = ^I'lSy V' • (4) 

Their corresponding charges, 

Qa=J d'x ^\xf-^^{x) , ^» ^ / '^'^ ^\x)'yj-^i^{x) , (5) 

are, likewise, generators of 5*17(2) x 5*17(2). 



2.2 Spontaneous chiral symmetry breaking 

There is evidence from hadron spectroscopy that the chiral SU{2) x 5(7(2) symme- 
try of the QCD Lagrangian Q with m = is spontaneously broken: for dynamical 
reasons of non-perturbative origin, the ground state (vacuum) of QCD has lost 
part of the symmetry of the Lagrangian. It is symmetric only under the subgroup 
SU{2)v generated by the vector charges . This is the well-known isospin sym- 
metry seen in spectroscopy and dynamics. 



Exploring the nucleus in the context of low-energy QCD 



5 



If the ground state of QCD were symmetric under chiral SU{2) x SU{2), both 
vector and axial charge operators would annihilate the vacuum: Q^|0) = Qa 10) = 0. 
This is the Wigner-Weyl realisation of chiral symmetry with a "trivial" vacuum. 
It would imply the systematic appearance of parity doublets in the hadron spec- 
trum. For example, correlation functions of vector and axial vector currents should 
be identical, i. e. {0| V'''V''|0) = {Q\A''A"\Q). Consequently, the spectra of vector 
(J'^ — 1~) and axial vector [J^ — 1^) mesonic excitations should also be identi- 
cal. This degeneracy is not seen in nature: the p meson mass (m^ ~ 0.77 GeV) is 
well separated from that of the a\ meson (rriai — 1.23 GeV). Likewise, the light 
pseudoscalar [J^ = 0~) mesons have masses much lower than the lightest scalar 
= 0+) mesons. 

One must conclude Qa\0} 7^ 0, that is, chiral symmetry is spontaneously broken 
down to isospin: SU{2)r x SU{2)l SU(2)v- This is the Nambu-Goldstone real- 
isation of chiral symmetry. Goldstone's theorem says that a spontaneously broken 
global symmetry implies the existence of a massless boson (the Goldstone boson). 
Let us give a sketch of this theorem. If Qa\Ci) 7^ 0, there must be a physical state 
generated by the axial charge, \(f>a) ~ Qa\0), which is energetically degenerate 
with the vacuum. Let Ho be the QCD Hamiltonian (with massless quarks) which 
commutes with the axial charge. Setting the ground state energy equal to zero 
for convenience, we have Ho\(t>a) ~ QaHo\0} = 0. Evidently \(j>a) represents three 
massless pseudoscalar bosons (for Nf = 2). They are identified with the pions. 

2.3 The chiral condensate 

Spontaneous chiral symmetry breaking goes together with a qualitative rearrange- 
ment of the vacuum, an entirely non-perturbative phenomenon. The ground state 
is now populated by scalar quark-antiquark pairs. The corresponding ground state 
expectation value (0|i/)'i/'|0) is called the chiral (or quark) condensate. We frequently 
use the notation 

iijij) = (uu) + (dd) . (6) 
The precise definition of the chiral condensate is: 

(tPtP) ^ -iTr lim SF{x,y) (7) 

with the full quark propagator, SF{x,y) = — i{0|Ti/'(a;)7/;(j/)|0) where T denotes 
the time-ordered product. We recall Wick's theorem which states that T^{x)il}{y) 
reduces to the normal product : ^{x)il}{y) : plus the contraction of the two field op- 
erators. When considering the perturbative quark propagator, Sp\x,y), the time- 
ordered product is taken with respect to a trivial vacuum for which the expectation 
value of : tfiil} : vanishes. Long-range, non-perturbative physics is then at the origin 
of a non-vanishing {: tp^p :). 

(In order to establish the connection between spontaneous chiral symmetry 
breaking and the non- vanishing chiral condensate in a more formal way, introduce 
the pseudoscalar operator Pa{x) = ip{x)'^5Taip{x) and derive the (equal-time) com- 
mutator relation [Qa,Pb] ~ —5ab'>pip which involves the axial charge Q'^ of eq.©. 
Taking the ground state expectation value, we see that Q^jO) 7^ is indeed consis- 
tent with {i!ip) / 0.) 
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Let |7ra(p)) be the state vectors of the Goldstone bosons associated with the 
spontaneous breakdown of chiral symmetry. Their four-momenta are denoted p** = 
{Ep,p), and we choose the standard normahzation 

{7Ta{p)\nt{p')} = 2Ep5at{27TfS^{p-p'). (8) 

Goldstone's theorem also implies non-vanishing matrix elements of the axial current 
(4) which connect \na{p)) with the vacuum: 

{0\A^,{x)\7v,{p)) = ip'*/5,,e-'''■^ (9) 

where / is the pion decay constant (taken here in the chiral limit, i. e. for vanishing 
quark mass). Its physical value 

U = (92.4 ± 0.3) MeV (10) 

differs from / by a small correction linear in the quark mass mq. 

Non-zero quark masses m^^d shift the mass of the Goldstone boson from zero 
to the observed value of the physical pion mass, m,r. The connection between m-n 
and the u- and d- quark masses is provided by PCAC and the Gell-Mann, Oakes, 
Renner (GOR) relation 

ml ^ ~^{mu+md){qq} +0{ml^a)- (H) 

We have set (qq) = (uu) ~ (dd) making use of isospin symmetry which is valid to 
a good approximation. Neglecting terms of order mJJ ^, identifying f = fn = 92.4 
MeV to this order and inserting rriu + rud — 12 MeV |17l I18| (at a renormalisation 
scale of order 1 GeV), one obtains 

{qq) ~ -(240 MeV)^ ~ -1.8 /m-^ (12) 

This condensate (or correspondingly, the pion decay constant /^) is a measure of 
spontaneous chiral symmetry breaking. The non-zero pion mass, on the other hand, 
reflects the explicit symmetry breaking by the small quark masses, with ~ mq. 
It is important to note that m, and {qq) are both scale dependent quantities. Only 
their product mq{qq) is invariant under the renormalization group. 

In order to appreciate the strength of the scalar condensate 11211 from a nuclear 
physics point of view, note that its magnitude is more than a factor of ten larger 
than the baryon density in the bulk of a heavy nucleus, p ~ 0.16 fm^"'. 

2.4 The nucleon mass and the gap in the hadron spectrum 

In the hadronic phase of QCD in which nuclei reside, confinement implies that the 
active ("effective") degrees of freedom are not elementary quarks and gluons but 
mesons and baryons. The observed spectrum of low-mass hadrons has a character- 
istic gap, A ~ Mjv ~ 1 GeV, which separates the masses of all baryons and almost 
all mesons from the ground state |0). On the other hand, the lightest pseudoscalar 
mesons are positioned well within this gap, for good reason: as Goldstone bosons 
of spontaneously broken chiral symmetry, pious would start out massless. Explicit 
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symmetry breaking by the masses of the light quarks introduces perturbations 
on a scale small compared to A. 

The appearance of the gap A is thought to be closely linked to the presence of 
the chiral condensate (tpip) in the QCD ground state. For example, loffe's formula 
|19|. based on QCD sum rules, connects the nucleon mass Mjv directly with {tpip} 
in leading order: 

MN = -^{^i>) + ... , (13) 

where Ab ~ 1 GeV is an auxiliary scale (the Borel mass) which separates "short" 
and "long" -distance physics in the QCD sum rule analysis. While this formula is 
not very accurate and needs to be improved by including higher order condensates, 
it nevertheless demonstrates that spontaneous chiral symmetry breaking plays an 
essential role in giving the nucleon its mass. 

The condensate (tptp) is encoded in the pion decay constant through the 
GOR relation IIH . In the chiral limit (niq — » 0), this is the only quantity which 
can serve to define a mass scale ("transmuted" through non-perturbative dynamics 
from the scale Aqcd — 0.2 GeV appearing in the QCD running coupling strength). 
It is common to introduce 47r/^ ~ 1 GeV as the scale characteristic of spontaneous 
chiral symmetry breaking. This scale is then roughly identified with the spectral 
gap A. 

At a more basic level, the nucleon mass is determined by the energy-momentum 
tensor T^" derived from the QCD Lagrangian Q . The component T"" = H defines 
the Hamiltonian, and Mjv = {N\H\N). A detailed analysis |20| shows that roughly 
2/3 of Mjv can be thought of as having its origin in the gluon field energy; the 
remainder results from the combination of quark kinetic and potential energies. 
Alternatively, one can express the nucleon mass as the expectation value Mjv = 
{N\Tl^\N) of the non- vanishing trace 

= ^G^.G"" -^m^uu + midd+ ... , (14) 

of the energy- momentum tensor (the so-called trace anomaly - see e.g. ref . [2T] ) . 
Here G*^" is the gluonic field strength tensor and fi[g) is the beta function of QCD 
(/? = —f3og^/ (47r)^ in leading order, with — 11 — 2Nf /3). The quark mass terms'^ 
mqqq with q = u, d . . . include light as well as heavier quarks. The physical nucleon 
mass A/jv is expressed as 

Mn = Mo + gn (15) 
in terms of its evidently non- vanishing value in the chiral limit, 

Mo = {N\-^G^,G^-' + ...\N) (16) 

2<? 

(where the dots refer to possible contributions from heavier quarks, other than u 
and d), and the sigma term 

crjv — mq^—^ = {N\muUU + mddd\N) . (17) 

q — u,d 



We have omitted here the anomalous dimension of the mass operator, as in |2()| . 
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This sigma term represents the contribution from exphcit chiral symmetry breaking 
to the nucleon mass, through the small but non-vanishing it- and d-quark masses. 

Relations such as Eg. HSU , when combined with Eg. HIGH , give important hints. 
Non-perturbative gluon dynamics confines quarks and generates hadron masses, 
producing a gap in the spectrum. In the sector of the (almost) massless u- and 
d- quarks, this gap reflects the spontaneously broken chiral symmetry. Systems 
characterized by an energy gap usually exhibit qualitative changes when exposed 
to variations of thermodynamic conditions. A typical example is the temperature 
dependence of the gap in a superconductor. For the physics in the hadronic phase 
of QCD, the following key issues need therefore to be addressed: how does the 
quark condensate (tp^p) change with temperature and/or baryon density? Are there 
systematic changes of hadronic spectral functions in a dense and hot medium, which 
would indicate changes of the QCD vacuum structure? How does the nucleon mass 
itself vary under such conditions? And what is the impact of these considerations 
on nuclear many-body systems? The sigma term 11711 will play an important role in 
this discussion as it controls the leading density dependence of the chiral condensate 
(see Section 2.7). 

2.5 Chiral effective field theory 

The mass scale given by the gap A ~ 47r/^ offers a natural separation between 
"light" and "heavy" (or, correspondingly, "fast" and "slow") degrees of freedom. 
The basic idea of an effective field theory is to introduce the active light particles 
as collective degrees of freedom, while the heavy particles are frozen and treated 
as (almost) static sources. The dynamics is described by an effective Lagrangian 
which incorporates all relevant symmetries of the underlying fundamental theory. 
We now list the necessary steps, first for the pure meson sector (baryon number 
B — 0) and later for the B = 1 sector. We work mostly with two quark fiavours 
{Nf = 2) in the following. 

a) The elementary quarks and gluons of QCD are replaced by Goldstone bosons. 
They are represented by a 2 x 2 matrix field U{x) £ SU{2) which collects the three 
isospin components TTaix) of the Goldstone pion. A convenient choice of coordinates 
is 

U{x) = exp[iTacl>a{x)] , (18) 
with (pa = T^a/f where the pion decay constant / in the chiral limit provides a 
suitable normalisation. (Other choices, such as U — -^/l — vTa//^ + iraiva/ f, are 
also common. Results obtained from the effective theory must be independent of 
the coordinates used.) 

b) Goldstone bosons interact weakly at low energy. In fact, if |7r) — (5^|0) is a 
massless state with H\n) = 0, then a state |7r") = (Q^)"|0) with n Goldstone bosons 
is also massless since the axial charges all commute with the full Hamiltonian H. 
Interactions between Goldstone bosons must therefore vanish at zero momentum 
and in the chiral limit. 

c) The QCD Lagrangian Q is replaced by an effective Lagrangian which in- 
volves the field U{x) and its derivatives: 



Cqcd c,ff{u,du,d^u,...). 



(19) 



Exploring the nucleus in the context of low-energy QCD 



9 



Goldstone bosons do not interact unless they carry non-zero four-momentum, so the 
low-energy expansion of 11911 is an ordering in powers of 9^(7. Lorentz invariance 
permits only even numbers of derivatives. We write 

= + (20) 

and omit an irrelevant constant term. The leading term, called non-linear sigma 
model, involves two derivatives: 

'Cl?'' = ^Trld^U'^d'^U]. (21) 

At fourth order, the terms permitted by symmetries are (apart from an extra con- 
tribution from the QCD anomaly, not included here): 

Z:^ = ^j.{Tr[d^U''d^'U]f + jTr[d^U''dM]Tr[d''U^d''U], (22) 

and so forth. The constants li, I2 (following canonical notations of ref. \'2'A}) must be 
determined by experiment. To the extent that the effective Lagrangian includes all 
terms dictated by the symmetries, the effective theory is equivalent to QCD |25j. 

d) The symmetry breaking mass term is small, so that it can be handled per- 
turbatively, together with the power series in momentum. The leading contribution 
introduces a term linear in the quark mass matrix m: 

£(.2) = L.Tr[d^u''d''U] + ^Bo Tr[m{U + (7^)]. (23) 

The fourth order term Ci^^ also receives symmetry breaking contributions with 
additional constants h. 

When expanding to terms quadratic in the pion field, one finds 

Cf^ = (m„ -I- ma)fBo + ad"" -k a - ^(m„ -f ma)Bo-Kl + 0(7r*). (24) 

At this point we can immediately verify the COR relation lllll at the level of the 
effective theory. The first (constant) term in 12411 corresponds to the shift of the 
vacuum energy density by the non-zero quark masses. Identifying this with the 
vacuum expectation value of the corresponding piece in the QCD Lagrangian Q, 
we find —{niuiuu) +md{dd)) = {mu + md)f^ Bo and therefore {uu) — {dd) — —f^Bo 
in the chiral limit, m^^d — > 0. The pion mass term in 1241 is evidently identified as 
= (m„ 4- md)Bo. Inserting Bq, we have the GOR relation. 

e) Given the effective Lagrangian, the framework for systematic perturbative 
calculations of the S-matrix involving Goldstone bosons, named Chiral Perturbation 
Theory (ChPT), is then defined by the following rules: 

Collect all Feynman diagrams generated by Ceff- Classify all terms according 
to powers of a variable Q which stands generically for three-momentum or energy 
of the Goldstone bosons, or for the pion mass m-n . The small expansion parameter 
is Q/Anf-n- Loops are subject to dimensional regularization and renormalization. 
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2.6 Effective Lagrangian including baryons 

The prominent role played by the pion as a Goldstone boson of spontaneously 
broken chiral symmetry has its strong impact on the low-energy structure and 
dynamics of nucleons as well. Decades of research in nuclear physics have established 
the pion as the mediator of the long-range force between nucleons. When probing 
the individual nucleon itself with long-wavelength electroweak fields, a substantial 
part of the response comes from the pion cloud, the "soft" surface of the nucleon. 

The calculational framework for this, baryon chiral perturbation theory, has 
been applied quite successfully to a variety of low-energy processes (such as thresh- 
old pion photo- and electroproduction and Compton scattering on the nucleon) 
for which increasingly accurate experimental data have become available in recent 
years. An introductory survey is given in chapter 7 of ref. |15|. A detailed review 
can be found in ref. |24l . 

Let us consider the sector with baryon number B — 1 and concentrate on the 
physics of the pion-nucleon system, restricting ourselves to the case of TV/ = 2 
flavours. 

The nucleon is represented by a Dirac spinor field <Fn{x) = {p,n)^ organised 
as an isospin-1/2 doublet of proton and neutron. The free field Lagrangian 

Co = {iif.d'' - Mo)<I'N (25) 

includes the nucleon mass in the chiral limit, Mq. Note that nucleons, unlike pions, 
supposedly have at least part of their large mass connected with the strong scalar 
mean field provided by the quark condensate {if'il'). Such a relationship is explicit, 
for example, in the loffe formula 11311 . 

We can now construct the low-energy effective Lagrangian for pions interacting 
with a nucleon. The previous pure meson Lagrangian is replaced by jC^f / = jCn + 
C^^^ + ... f which also includes the nucleon field. The additional term involving 
the nucleon is expanded again in powers of derivatives (external momenta) and 
quark masses: 

rf„ = C^l + C^:i ... (26) 

Let us discuss the leading term, -C^^^. The modifications as compared to the free 
nucleon Lagrangian (I25I I are twofold. First, there is a replacement of the d'^ term by 
a chiral covariant derivative which introduces vector current couplings between the 
pions and the nucleon. Secondly, there is an axial vector coupling. This structure 
of the TvN effective Lagrangian is again dictated by chiral symmetry. We have 

Ci^l ^ ^Nlij^d^ - iv'') + -y^-y^a'' - MoWn , (27) 

with vector and axial vector quantities involving the Goldstone boson (pion) fields 
in the form ^ = \/TJ : 

^ l(^t£,M^ ^ ^Q^^^t■^ ^ a'Vc + ... , (28) 

«" = ^(e^a''e-C9"e^) = -^r„a''^„-f... , (29) 

where the last steps result when expanding v'^ and a'' to leading order in the pion 
fields. 
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So far, the only parameters that enter are the nucleon mass, Mo, and the pion 
decay constant, /, both taken in the chiral limit and ultimately connected with a 
single scale characteristic of non-perturbative QCD and spontaneous chiral sym- 
metry breaking. 

When adding electroweak interactions to this scheme, one observes an additional 
feature which has its microscopic origin in the substructure of the nucleon, not 
resolved at the level of the low-energy effective theory. The analysis of neutron 
beta decay (n — > pev) reveals that the 7/^75 term in 12711 is to be multiplied by the 
axial vector coupling constant qa, with the empirical value 

gA = 1.267 ± 0.003 . (30) 

At next-to-leading order {C''^^), the symmetry breaking quark mass term enters. 
It has the effect of shifting the nucleon mass from its value in the chiral limit to 
the physical one, introducing the sigma term ajv of Eg. 11711 . The empirical sigma 
term, 

(jjv = (45 ± 8)MeV , (31) 

has been deduced by a sophisticated extrapolation of low-energy pion-nucleon data 
using dispersion relation techniques 1251 . There is still an ongoing debate, related 
to different ways of handling the ttN phase shift analysis, whether crjv could in fact 
be larger than that given in 13H . On the other hand, a recent extrapolation of 
two-flavour lattice QCD results using methods of chiral effective field theory |26| 
yields" aN = (47 ± 3) MeV, perfectly consistent with Eg.CT. 

Up to this point, the nN effective Lagrangian, expanded to second order in the 
pion field, has the form 

C^ff = >^iv(j7M9^ - Mn)<Pn - ^^N-f^ToT^N ■ a^TT (32) 

2/n- 

4/77 /tt 

where we have not shown a series of additional terms of order (d'^Tv)^ included in 
the complete jC]^j^- These terms come with further constants that need to be fitted 
to experimental data. In later applications using in-medium chiral perturbation 
theory, the calculations will be based on the leading terms, linear in the derivative 
d^n of the pion field, while the term involving crjv will be manifest in the density 
dependence of the chiral condensate, as we shall see now. 



2.7 Chiral thermodynamics 

Before turning our attention to the nuclear many-body system, it is of some rel- 
evance to examine the thermodynamics based on the chiral effective Lagrangian 
^eff = ~^ ^ttn)- particular, we are interested in the question how the 

quark condensate (qq) behaves with changing thermodynamic conditions. 

* note: statistical errors only; systematic lattice uncertainties from finite volume 
effects etc. still require further detailed studies. 
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A. Thermodynamics of the chiral condensate 

Consider as a starting point the partition function 

1 



Z — Tr exp 



(33) 



where Tl is the Hamiltonian density, fi denotes the chemical potential, and p the 
baryon density. The partition function a,t p — 0, expressed in terms of the QCD 
Hamiltonian H, is 

Z = Trexp{-H/T) = ^{n\e''^"^'^\n) (34) 

n 

with {H — En)\n) = 0. Confinement implies that the eigenstates |n) of _ff are 
(colour-singlet) hadrons at low temperatures, below the critical Tc ~ 0.2 GeV for 
the deconfinement transition. The physics is then determined by the states of low- 
est mass in the spectrum {En}, namely the pions. At non-zero baryon chemical 
potential p, nucleons enter in addition as the lowest-mass baryons. The Hamilto- 
nian density Tl of QCD is expressed in terms of the relevant low-energy degrees of 
freedom in the hadronic phase, derived from the chiral effective Lagrangian C^ff- 
The quark mass term, 

STi. — tprmf} — rriu uu + TUd dd + ... , (35) 

can be separated so that Ti. = Tio + STl, with Tio representing the massless limit. 

Now assume a homogeneous medium and consider the pressure (or the free 
energy density) 

P{T,V,p) = ~HT,V,p) = ^lnZ . (36) 

The derivative of P with respect to a quark mass of given flavour q = u,d 
obviously produces the in-medium quark condensate, the thermal expectation value 
{qqjr as a function of temperature and chemical potentials. Subtracting vacuum 
quantities one finds 

/- \ /- \ dP{T,V,p) 

{qq)T,p = (gg)o J , (37) 

where (gg)o refers to the vacuum condensate taken at T = and ^ = 0. The 
/i-dependence of the condensate is converted into a density dependence via the 
relation p = dP/dp at fixed T. Using the Cell-Mann, Oakes, Renner relation lllll . 
one can rewrite Eg. 13711 as 

{gg)T.p _ 1 I 1 dP{T,p) ^ 
{ot)o /I dml 

The task is therefore to investigate how the equation of state changes, at given tem- 
perature and baryon chemical potential, when varying the quark mass (or, equiva- 
lently, the squared pion mass). 

For a system of nucleons interacting with pions, the total derivative of the 
pressure with respect to m\ reduces to 

dP dP omm dP dP OM . 



dmi dmi dm% OMn drn^ ' 
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with the scalar density ps = —dP/dMN and the sigma term crjv = m%dMN / drn^ 
translated from Eg. 1)17^ using the GOR relation. 

Next, consider the limit of low baryon density p at zero temperature, T = 0. 
In this case nucleons are the only relevant degrees of freedom. At sufficiently low 
density, i.e. when the average distance between two nucleons is much larger than the 
pion Compton wavelength, the nucleon mass remains at its vacuum value A/jv , and 
one can neglect A'^A'^ interactions. The pressure is that of a Fermi gas of nucleons^, 
subject only to the Pauli principle. Returning to Eas. l38l39ll . we have 

{qq)o mif4 

with the scalar density 



P=-^=a[ '^'p (41) 



for a system with equal number of protons and neutrons, i.e. with degeneracy factor 
d = 4 from spin and isospin. The Fermi momentum kf is related to the baryon 
density by p = 2fc3/(37r^). At low baryon densities with <^ M%, the difference 
between ps and p can be neglected, so that 

{qq)o mlfi 

The consequences implied by this relation are quite remarkable. Using the empir- 
ical value aN — 50 MeV of the nucleon sigma term one observes that the chiral 
condensate at normal nuclear matter density, p = po = 0.16 fm~^, is expected to 
decrease in magnitude to less than 2/3 of its vacuum value, a significant effect that 
should have observable consequences already in the bulk parts of ordinary nuclei. 

In the limit of low temperature and at vanishing baryon density (p = 0), pious 
are the only thermally active degrees of freedom in the partition function H34|l . This 
case is well described by chiral perturbation theory. The combined low-temperature 
and low-density behaviour is summarised as 

(W)p _ 1 Q 3^ / crjv \ p 



{qq)o 8/2 • \50MeVj po 

showing a rather weak leading dependence of the chiral condensate on temperature, 
whereas its density dependence is far more pronounced. 



B. A schematic model 

Consider now a schematic model for the hadronic phase of QCD, starting from an 
effective Lagrangian Ceff ~ Ctt + iCttN + jCnn in which £^ is given by Ea. ll24ll . C^n 
is the pion- nucleon Lagrangian of Eg. 1)27^ and additional short distance interactions 
between nucleons are represented by A'^A*' contact terms, 

^ Actually, this statement needs to be modified at extremely low densities where 
the system turns into a gas of clusters (deuterons etc.) 
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Fig. 1. Chiral condensate (in units of its vacuum value) as a function of temperature 
and baryon density (po = 0.16 fm""^ is the density of normal nuclear matter). 



Cnn = ^{^N<pNf - ^('^ATTM^'iv)^ + - ■ (44) 

Their coupling strength parameters Gs < and Gy > are fixed to reproduce 
ground state properties of normal nuclear matter. What we have in mind here is 
a variant of relativistic mean field theory combined with "soft" pion fluctuations 
treated in the framework of chiral perturbation theory and to be discussed in more 
detail later. 

Using two-loop thermal field theory in order to perform a self-consistent calcu- 
lation of the pressure P{T, fi) in this model, one can deduce the chiral condensate as 
a function of temperature and baryon density following Eg. 1381 1. This calculation 
generates temperature dependent mean fields for the nucleons at the same time 
as it treats thermal pion fiuctuations with leading nn interactions. The pressure 
equation takes the form 

P{T,fi)=PN{T,fi*,M*^)+P^{T,,,*,Ar^) + ^pi + ^p^ , (45) 

with nucleon and pion contributions Pn,tv respectively depending on the effective 
nucleon mass M'^ = Mn + Gs ps and the shifted baryon chemical potential p* = 
p — Gv p. The baryon and scalar densities are determined as 

_dP _ djPN + P^) _ dP _ djPN + Pn) 

dp' dp* ' P'- omn ' dM*^ ■ ^^^> 

Fixing Gs,v to the energy per particle and the equilibrium density of cold nuclear 
matter, one can map out the nuclear equation of state first at low temperature 
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and density, reproducing known physics in this domain which is a prerequisite for 
extrapolating into more extreme regions. 

We now return to the evaluation of the chiral condensate. The dependence 
of P{T, n) on the pion mass is explicit in the thermal pion Green function and 
implicit through the nucleon mass. The result |27| for {qq)T,p is shown in Fig0 At 
low baryon density the linear behaviour of Eg. 1421 is recovered. Pionic fluctuations, 
calculated up to three-loop order with inclusion of two-pion exchange effects, help 
maintain this approximately linear dependence not only up to p ~ po, but even 
slightly beyond. Up to this point we can conclude that the magnitude of the quark 
condensate at normal nuclear matter density is expected to be reduced by about 
one third from its vacuum value, whereas the temperature dependence is far less 
pronounced, at least up to T < 100 MeV. 



3 Chiral dynamics and the nuclear many body problem 

Nuclei are aggregates of nucleons (and mesons) which are in turn clusters of quarks 
and gluons. Nuclei are thus central cornerstones in the QCD phase diagram. The 
decription of nuclear many-body systems must ultimately be linked to and con- 
strained by QCD. Recent developments using effective field theory methods are 
aiming to establish such connections. We focus here on a framework that relates 
strongly to principles of spontaneous chiral symmetry breaking in QCD and to the 
topics presented in the previous sections. 

We propose a two-step approach. The first step emphasises the role of pions 
in nuclear many-body dynamics, given the rules dictated by chiral symmetry. Our 
strategy is based on in-medium chiral perturbation theory where the nuclear Fermi 
momentum enters as an additional scale in the systematic low-energy expansion 
of chiral effective field theory. As will be demonstrated, this approach is already 
capable of producing binding and saturation of nuclear matter. At the same time 
it gives a remarkably good value for the asymmetry energy with no additional 
parameter. At that level, however, important features such as the nuclear spin- 
orbit force cannot yet be understood. The second step is to provide strong scalar 
and vector mean fields, about equal in magnitude but opposite in sign, which are 
known phenomenologically to account for the large spin-orbit splitting in finite 
nuclei. These strong mean fields will be viewed as arising from the changes of the 
condensate structure of the QCD vacuum in the presence of the nuclear medium. 
QCD sum rules at finite baryon density are used for guidance to set constraints on 
the in-medium modifications of the relevant condensates. 

3.1 Nuclear matter, part I 

The present status of the nuclear matter problem can briefiy be summarised as fol- 
lows: a quantitatively successful description is achieved, using advanced many-body 
techniques |28|. in a non-relativistic framework when invoking an adjustable three- 
body force. Complementary to such an ab-initio approach, relativistic mean field 
phenomenology, including non-linear terms with adjustable parameters, have also 
been widely used for the calculation of nuclear matter properties and finite nuclei 0. 
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At a more basic level, the Dirac-Brueckner method |29| solves a relativistically im- 
proved Bethe-Goldstone equation with one-boson exchange A'^A'^-interactions. The 
resulting Brueckner G-matrix is then used as input in Hartree-Fock calculations. 

In recent years a novel approach to the NN interaction based on effective field 
theory (in particular, chiral perturbation theory) has emerged 1301 1311 1^ I33| . Its 
key element is a power counting scheme which separates long- and short-distance 
dynamics. Methods of effective field theory have also been applied to systems of 
finite density ^Hl^ . 

In the following sections the primary purpose is to point out the importance 
of explicit pion dynamics in the nuclear many-body problem. While pion exchange 
processes are well established as generators of the long and intermediate range 
A'^A'" interaction, their role in nuclear matter is less evident. The one-pion exchange 
Hartree term vanishes identically for a spin-saturated system, and the leading Fock 
exchange term is small. Two-pion exchange mechanisms are commonly hidden be- 
hind a purely phenomenological scalar ("sigma") field which is fitted to empirical 
data but has no deeper justification. This situation calls for a more detailed un- 
derstanding. We report on steps in this direction, following ref. 9 • A conceptually 
similar strategy has been proposed in ref. 36 . 

Before passing on to explicit calculations it is useful to draw attention to the 
following fact. A simple but surprisingly good parametrization of the energy per 
particle, E{kf) — E{kf)/A, of isospin symmetric nuclear matter is given in powers 
of the Fermi momentum fc/ as 

where the nucleon density is p = 2fc^ /37r^ as usual, and M jv is the free nucleon mass. 
The first term is the kinetic energy of a Fermi gas. Adjusting the (dimensionless) 
parameters a and /3 to the equilibrium density, po = 0.16 fm""^, (fc/o = 1.33 fm~^) 
and _Eo — E{kfQ) = — 16 MeV, gives a = 5.27 and (3 = 12.22. The compression 
modulus Ko — k^Q{d^ E{kf)/dk^)kjQ is then predicted at Kq = 236 MeV, well in 
line with empirically deduced values, and the density dependence of E{kf) using eq. 
14711 is remarkably close to the one resulting from the realistic many-body calcula- 
tions of the Urbana group |37| . The expansion 1471 in powers of Fermi momentum 
(rather than in powers of density) is the natural one in the context of in-medium 
chiral effective field theory to which we now turn our attention. 

3.2 In-medium chiral perturbation theory 

The tool to investigate the implications of spontaneous and explicit chiral symmetry 
breaking in QCD is chiral perturbation theory. Observables are calculated within 
the framework of an effective field theory of Goldstone bosons (pions) interacting 
with the lowest-mass baryons (nucleons), as described in sections 2.5 and 2.6. The 
diagrammatic expansion of this low-energy theory in the number of loops has a one- 
to-one correspondence to a systematic expansion of observables in small external 
momenta and the pion (or quark) mass. 

In nuclear matter, the relevant momentum scale is the Fermi momentum fc/. 
At the empirical saturation point, fc/o — 2m^, so the Fermi momentum and the 
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pion mass are of comparable magnitude at the densities of interest. This imme- 
diately implies that pions should be included as explicit degrees of freedom: their 
propagation in matter is relevant. Pionic effects cannot be accounted for simply by 
adjusting coefficients of local A''A'' contact interactions. 

Both kf and are small compared to the characteristic chiral scale, 47r/^ ~ 
1.2 GeV. Consequently, the equation of state of nuclear matter as given by chiral 
perturbation theory will be represented as an expansion in powers of the Fermi 
momentum. The expansion coefficients are non-trivial functions of kf/m-n, the di- 
mensionless ratio of the two relevant scales inherent to the problem. 

The chiral effective Lagrangian 1271 generates the basic pion-nucleon coupling 
terms: the Tomozawa- Weinberg tvttNN contact vertex, (l/4/^)(g^ — qa)'yfitabcTc, 
and the pseudovector ttNN vertex, (5A/2/,r)'7a7M75''"a, where qa,b denote (outgoing) 
pion four- momenta and qa is the axial vector coupling constant (we choose qa ~ IS 
so that the Goldberger-Treiman relation Qt^n = jtaM//^ gives the empirical nN 
coupling constant, Qt^n = 13.2). The only new ingredient in performing calculations 
at finite density (as compared to evaluations of scattering processes in vacuum) is 
the in-medium nucleon propagator. For a relativistic nucleon with four-momentum 
p** — {po,p) it reads 



The second term is the medium insertion which accounts for the fact that the 
ground state of the system has changed from an "empty" vacuum to a filled Fermi 
sea of nucleons. Diagrams can then be organized systematically in the number of 
medium insertions, and an expansion is performed in leading inverse powers of the 
nucleon mass, consistently with the A:/-expansion. 

Our "inward-bound" strategy is now as follows. One starts at large distances 
(small kf) and systematically generates the pion-induced correlations between nu- 
cleons as they develop with decreasing distance (increasing kf). The calculations 
of the energy density reported here are performed to 3-loop order (including terms 
up to order fcj) and incorporate one- and two-pion exchange processes. The pro- 
cedure involves one single momentum space cutoff A which encodes dynamics at 
short distances not resolved explicitly in the effective low-energy theory. This high- 
momentum scale A is the only free parameter at this stage. It has to be fine-tuned. 
(Alternatively and equivalently, one could use dimensional regularization, remove 
divergent loop integrals and replace them by adjustable A^'A^ contact terms which 
then parametrize the unresolved short-distance physics). 

Nuclear chiral dynamics up to three-loop order introduces the diagrams, Fig.|5| 
They include the one-pion exchange (OPE) Fock term, iterated OPE and irreducible 
two-pion exchange. Medium insertions are systematically applied on all nucleon 
propagators, and the relevant loop integrations yield results which can be written 
in analytic form for all pieces. 

We now outline the leading contributions to the energy per particle E(kf). The 
kinetic energy including first order relativistic corrections is 




(48) 




(49) 
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Fig. 2. In-medium chiral perturbation theory: One-pion exchange Fock term (top 
left), iterated one-pion exchange (top middle and right) and examples of irreducible 
two-pion exchange terms (bottom). See Ref. for details. 



Terms of order are already negligibly small. At least from this perspective, 
nuclear matter is a non-relativistic system. 

The OPE Fock term (see the diagram on the left of Fig. |2l becomes 



(47r/.)2 



F 



fc/ ^ , t k 



m 



(-) 



(50) 



where F and G are functions of the dimensionless variable fc//m^. They are given 
explicitly in Ref. W. All finite parts of iterated OPE and irreducible two-pion 
exchange are of the generic form 



(47r/,r) L Vm^/ Vm^ 



(51) 



with the functions -ff4,5 again given explicitly in Ref.jH]- AH power divergences 
specific to cutoff regularization are summarized in the expression 



EA{kf) = 



Akj 



(52) 



where the attractive and dominant first term in the brackets arises from iterated 
OPE. Note that this term, proportional to the density p ^ kj, can be generated in 
an equivalent mean-field approach by a A'^A'^ contact interaction with appropriate 
coupling strength. The sum of the terms H49^ - I|52|l determines the energy per particle 
of nuclear matter within chiral dynamics, to the order shown in Fig. |21 



3.3 Nuclear matter equation of state 

A striking feature of the chiral dynamics approach is the simplicity of the saturation 
mechanism for isospin-symmetric nuclear matter. Before turning to the presentation 
of detailed results, it is instructive first to discuss the situation in the exact chiral 
limit, ra-n = 0. The basic saturation mechanism can already be demonstrated by 
truncating the one- and two-pion exchange diagrams at order We can make 
straightforward contact with the parametrization Ij47|l of the energy per particle 
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Fig. 3. Energy per particle, E{kf), of symmetric nuclear matter derived from chiral 
one- and two-pion exchange (solid line) j^. The cutoff scale is vl = 646 Me V^. The 
dashed line is the result of Ref. |H7) . 

and identify the coefficients a and /3 of the kj and fc* terms, respectively. The 
result for a in the chiral limit is: 




where we have neglected the small correction proportional to in Eg. 152^ . The 
strongly attractive leading term in Ea. H52^ is accompanied by the (weakly repulsive) 
one-pion exchange Fock term. The fc|-contribution to E{kf) would lead to collapse 
of the many-body system. The stabilizing fc^-term is controlled by the coefficient 
(calculated again in the chiral limit) 

^ = 4 (^)' + 237 - 24 In 2) - I = 13.55, (54) 

a unique and parameterfree result to this order. The two-pion exchange dynamics 
in combination with the Pauli exclusion principle produces repulsion of just the 
right magnitude to achieve saturation: the result, Ea. (l54L is within 10% of the 
empirical (3 — 12.2. Adjustment of the short-distance scale A between 0.5 and 0.6 
GeV easily leads to a stable minimum of E{kf) in the proper range of density and 
binding energy. 

The full 3-loop chiral dynamics result for E{k;) in symmetric nuclear matter, 
using = 135 MeV (the neutral pion mass), is shown in Fig|3| together with a 
realistic many-body calculation |37|. The outcome is remarkable: with one single 
parameter A = 0.65 GeV fixed to the value Eo — —15.3 MeV at equilibrium. 
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perturbative pion dynamics alone produces an equation of state which follows that 
of much more sophisticated calculations up to about three times the density of 
nuclear matter. The predicted compression modulus is Kq = 255 MeV, well in line 
with the "empirical" Kq — (250 ± 25) MeV. The basic mechanisms behind binding 
and saturation in this approach are thus identified as follows. The exchange of two 
pions between nucleons (truncated at the momentum scale A) drives the attraction'^ 
proportional to fcj in the energy per particle. The Pauli blocking of intermediate 
nucleon states in these 27r exchange processes produces the (repulsive) stabilizing 
k'f term in E/A. 

The calculations just described have also been extended to finite temperature 
1^. The resulting pressure as function of density shows the van der Waals - like 
behaviour of the liquid - gas phase transition, with a critical temperature Tc — 
25 MeV, slightly larger than the value around 18 MeV commonly quoted as the 
"empirical" one. 

3.4 Asymmetry energy 

The specific isospin dependence of two-pion exchange should have its distinct 
influence on the behaviour of asymmetric nuclear matter, with increasing ex- 
cess of neutrons over protons. We introduce as usual the asymmetry parameter 
S = {pn — Pp)/p = {N — Z)/{N+Z), keeping the total density p = p„+Pp — 2kf /3tv^ 
constant. The proton and neutron densities are pp^„ = {k^'")^ /3n^ in terms of the 
corresponding Fermi momenta, fc^ and kf. Without change of any input, the asym- 
metry energy A{kf) defined by 

Easik"}, fc^) = E{kf) + S^A{kf) + ... (55) 

has been calculated 0. The result at nuclear matter density is Ao = A{kfo) — 33.8 
MeV. This is in very good agreement with the empirical value Ao — (30 ± 4) MeV 
derived from extensive fits to nuclide masses |38| . 

Extrapolations to higher density work roughly up to p ~ 1.5 po- At still higher 
densities, there are indications that non-trivial isospin dependence at shorter dis- 
tances, beyond one- and two-pion exchange, starts to play a role. A similar state- 
ment holds for pure neutron matter which comes out properly unbound, but its 
predicted equation of state starts to deviate from that of realistic many-body cal- 
culations at neutron densities larger than 0.2 fm~^. 

3.5 Nuclear mean field from chiral dynamics 

The in-medium three-loop calculation of the energy per particle defines the (mo- 
mentum dependent) self-energy of a single nucleon in nuclear matter up to two-loop 
order. The real part of the resulting single particle potential in isospin-symmetric 
matter at the saturation point, for a nucleon with zero momentum, comes out as 

nsi 

U{p = 0, kfo) = -53.2 MeV, (56) 

® In a qualitative (but not strict) analogy, these mechanisms are reminiscent of 
van der Waals forces. 
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using exactly the same one- and two-pion exchange input that has led to the solid 
curve in Fig.|3| The momentum dependence of U{p, fc/o) can be rephrased in terms 
of an average effective nucleon mass M* ~ 0.8 Mn at nuclear matter density, and 
the imaginary part of the potential for a nucleon-hole at the bottom of the Fermi 
sea is predicted to be about 30 MeV. While these numbers are remarkably close to 
the empirically deduced ones, it must nevertheless be emphasised that the detailed 
momentum dependence of the single particle potential, as obtained in this approach, 
is not yet satisfactory. 

3.6 Intermediate summary and a working hypothesis 

Explicit pion dynamics originating from the spontaneously broken chiral symme- 
try of QCD is clearly an important aspect of the nuclear many-body problem. 
In-medium chiral perturbation theory, with one single cutoff scale A ~ 0.65 GeV 
introduced to regularize the few divergent parts associated with two-pion exchange, 
gives realistic binding and saturation of nuclear matter already at three-loop or- 
der. At the same time it gives very good values for the compression modulus and 
the asymmetry energy. These are non-trivial observations, considering that it all 
works with only one adjustable parameter encoding short-distance dynamics which 
remains unresolved at the Fermi momenta kf <^ Mn ~ 47r/,r characteristic of 
equilibrium nuclear matter. 

In view of the relevant scales in nuclear matter, the importance of explicit 
pion degrees of freedom does not come unexpected. Many of the existing models 
ignore pious, however. They introduce purely phenomenological scalar fields with 
non-linear couplings and freely adjustable parameters in order to simulate two-pion 
exchange effects. 

Now, in order to establish contacts with the widely used and successful relativis- 
tic nuclear mean field phenomenology, the following working hypothesis suggests 
itself as a guide for the next necessary step. Assume that the nuclear matter ground 
state represents a "shifted" QCD vacuum characterized by strong vector (V) and 
scalar (S) condensate fields acting on the nucleons, with V ~ —S ~ 0.3 GeV as 
suggested e. g. by in- medium QCD sum rules 013 |H|. Such a scenario would not 
produce binding all by itself. Binding and saturation would instead result from pio- 
nic (chiral) fluctuations as outlined previously. The strong condensate fields would, 
on the other hand, generate the large spin-orbit splitting in finite nuclei. This is 
the conceptual framework that we now apply when turning to finite nuclei. 

4 Relativistic nuclear model constrained by QCD and 
chiral symmetry 

At this point it is appropriate to recall in greater detail the relativistic nuclear phe- 
nomenology that goes under the name of Quantum Hadrodynamics (QHD) [J: a 
framework of Lorentz-covariant meson-baryon effective field theories. In the mean- 
field approximation, such an approach is equivalent to a model with local four-point 
interactions between nucleons Models based on QHD have been suc- 

cessfully applied to describe a broad range of nuclear phenomena, from light nuclei 
to superheavy elements (see Ref. '2' for a recent review, and references therein). 
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This framework has also been extended to studies of the structure of exotic nuclei 
with extreme isospin and close to the particle drip lines. 

The effective hadronic Lagrangians of QHD and related models incorporate 
known long-range interactions constrained by symmetries and a set of generic 
short-range interactions. These models are consistent with fundamental symmetries: 
Lorentz invariance, parity invariance, electromagnetic gauge invariance, isospin and 
with the chiral symmetry of QCD [Tll51l42|. The most successful QHD models are, 
however, purely phenomenological, with parameters adjusted to reproduce the nu- 
clear matter equation of state and global properties of spherical closed-shell nuclei. 
QHD calculations do not include pions explicitly, whereas we have just demon- 
strated that pions, as Goldstone bosons of spontaneously broken chiral symmetry, 
play an important role in the nuclear many-body problem. Two-pion exchange 
effects are supposedly incorporated as part of the strong scalar-isoscalar field of 
QHD models, but in an ad-hoc manner without detailed reference to the underly- 
ing nnNN dynamics. QCD symmetries do constrain the effective QHD Lagrangians 
by restricting the form of possible interaction terms. However, the empirical data 
set of bulk and single-particle properties of finite nuclei can only determine six or 
seven parameters in the general expansion of the effective Lagrangian in powers of 
the fields and their derivatives (43 . Such a general expansion can be controlled by 
the "naive dimensional analysis" (NDA) NDA tests the coefficients 

of the expansion for "naturalness" , i.e. this method controls the orders of magnitude 
of the coupling constants. NDA can exclude some interaction terms because their 
couplings would be "unnatural" , but it cannot determine the model parameters at 
the level of accuracy required for a quantitative analysis of nuclear structure data. 

Our approach is similar in spirit but proceeds with a different strategy, imposing 
as many low-energy QCD constraints as possible in order to minimize the number 
of freely adjustable parameters, ft is based on the conjectures already presented in 
the Introduction, namely: 

that nuclear binding and saturation arise dominantly from chiral (pionic) fluc- 
tuations in combination with Pauli blocking effects, as described in the previous 
section; 

that these pionic fluctuations are superimposed on individually strong scalar 
and vector fields of approximately equal magnitude and opposite sign, originating 
from density-dependent changes of the QCD vacuum condensates. 

Our aim is thus to study the interplay between condensate background fields 
and perturbative chiral fluctuations, both rooted in the spontaneous symmetry 
breaking pattern of QCD, in forming nuclei. We will demonstrate how this scenario 
works at large. Whereas in first approximation the condensate fields do not play a 
significant role for the saturation mechanism, they are essential for the description 
of ground states of finite nuclei and, in particular, the spin-orbit force. 

4.1 Point-coupling model: density dependent contact interactions 

A suitable framework to address these questions, both for nuclear matter and for 
finite nuclei, is a relativistic point coupling model with density dependent A'^A'' 
contact interactions. The explicit density dependence of these couplings reflects 
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effects "beyond mean field" such as they arise from the Pauli corrections to two-pion 

exchange fluctuations discussed previously. The mapping of the kf expansions of 
chiral dynamics, described in Section 3, into the form of equivalent point couplings 
depending on fractional powers of p, is performed by equating the nucleon self- 
energies resulting from both schemes. 

A. Lagrangian 

The relativistic point-coupling Lagrangian is built from basic density and current 
operator blocks bilinear in the Dirac spinor field ^ of the nucleon'^ : 

i^Orm , Or€{l,Ti} , r€{lni^,75,75Jl^,<ri^^}- (57) 

Here n are the isospin matrices and F generically denotes the Dirac matrices. The 

interaction terms of the Lagrangian are products of these blocks to a given order. In 
principle, a general effective Lagrangian can be written as a power series in the point 
couplings {'J^OtTW) and their derivatives. In practice, properties of symmetric and 
asymmetric nuclear matter, as well as empirical ground state properties of finite 
nuclei, constrain only the isoscalar-scalar, the isoscalar-vector, the isovector-vector, 
and to a certain extent the isovector-scalar channels. 

The model that we consider here includes the following four-fermion interaction 
vertices: 

isoscalar-scalar: {WWY 
isoscalar-vector: {9-f^ip){^')'^9) 
isovector-scalar: {^t9) ■ {'Pt^) 
isovector-vector: (#T7^*') • (#r7^*') . 

Vectors in isospin space are denoted by boldface symbols. The model is defined by 

the Lagrangian density 

•C = Cficc + + -Cder + -Cam, (58) 

with the four terms specified as follows: 

Ue. = ^{ij^d^ - M)9 , (59) 
= -i Gs(p)(#>f)(#>f) 

-\ Gts{pWt^) ■ {^T'^) 

- \ Gtv (p) (!^r7^<f ) • (<^r7''<f ) , (60) 

Cicr = -i Ds(/5)(a#*')(r#f) 

^ From here on we denote the nucleon field by dropping the index A'^. Similarly, 
the nucleon mass is denoted M without further index specification. 
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"2 ^Tv(/5)(9,.#T7^*') • (r#T7^<p-) , 



1 
2 
1 
2 
1 
2 



I>v'(/3)(9.#7m*')(o'"''^7''*') 



(61) 



(62) 



While the Lagrangian 1581 - 162(1 is going to be formally used in the mean-field ap- 
proximation, fluctuations beyond mean field are nevertheless encoded in the density- 
dependent couplings Gi{p) and Di{p), to be specified in detail later. When applied 
to finite nuclei the model must also include the coupling Lcm of the protons to 
the electromagnetic field, and the derivative terms in yCder, in addition to the free 
nucleon Lagrangian jCfree and the interaction terms in £45 . The derivative terms are 
built from the blocks dvi^^TiFj^l/) (the derivative is understood to act on both # 
and 'I'). One can, of course, construct many more derivative terms of higher orders. 
In £dcr we include only those terms which are direct counterparts of the second or- 
der interaction terms incorporated in Cm- In this way we take into account leading 
effects of finite range interactions that are important for a quantitative fine-tuning 
of nuclear properties. 

Our model Lagrangian formally resembles the ones used in the standard rela- 
tivistic mean-field point-coupling models of Refs. [89114011111 145II46| . The underlying 
dynamics is, however, quite different. The parameters of the interaction terms in 
the standard point-coupling models have constant values adjusted to reproduce the 
nuclear matter equation of state and a set of global properties of spherical closed- 
shell nuclei. In order to describe properties of finite nuclei on a quantitative level, 
these models include also some higher order interaction terms, such as six-fermion 
vertices (#<Z')^, and eight-fermion vertices (#*')'' and [(>^7^<Z')(>^7''<Z')]^. Our ap- 
proach includes only second order interaction terms, but with coupling strengths 
that are functions of the nucleon density operator p. Their functional dependence 
will be determined by matching the nucleon self-energies resulting from in-medium 
chiral perturbation theory and finite-density QCD sum rules. 

Medium dependent vertex functions have also been considered in so-called 
density dependent relativistic hadron field (DDRH) models. The density depen- 
dent meson-nucleon vertex functions are determined either by mapping the nu- 
clear matter Dirac-Brueckner nucleon self-energies in the local density approxima- 
tion 1471 1481 1491 . or they are adjusted to reproduce properties of symmetric and 
asymmetric nuclear matter and finite nuclei |50l I51| . In practical applications of 
the DDRH models the meson-nucleon couplings are assumed to be functions of the 
baryon density In a relativistic framework the couplings can also depend on 

the scalar density '1'^. Nevertheless, expanding in •f'^'f' is the natural choice, for 
several reasons. The baryon density is connected to the conserved baryon number, 
unlike the scalar density for which no conservation law exists. The scalar density 
is a dynamical quantity, to be determined self-consistently by the equations of mo- 
tion, and expandable in powers of the Fermi momentum. For the meson exchange 
models it has been shown that the dependence on baryon density alone provides a 
more direct relation between the self-energies of the density- dependent hadron field 
theory and the Dirac-Brueckner microscopic self-energies |49| . Moreover, the pion- 
exchange contributions to the nucleon self-energy, as calculated using in-medium 
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chiral perturbation theory, are directly given as expansions in powers of the Fermi 
momentum. Following these considerations, we express the coupling strengths of the 
interaction terms as functions of the baryon density, represented by the operator 
p = 1?/+!// in the rest frame of the many-body system. 



B. Equation of motion and nucleon self-energies 

The single-nucleon Dirac equation is derived by variation of the Lagrangian I58I I 
with respect to W: 

[7^(ir-V'^)-(A/ + S)]*- = 0, (63) 

where 

= + r ■ X'^ + + r ■ S'^^j, (64) 

and 

S = Es+T- StS + ErS + T ■ SrTS , (65) 

with the nucleon self-energies defined by the following relations 

= (Gv - DvO)f ~ eyl"-!-^ (66) 

S'i, ^ {Gtv - DTvO)j''^ (67) 
Ss ^ {Gs - Dsa){^1') (68) 
Sts = (Gts - DtsO){^t<]/) (69) 

Srs = -^{d.fWid'^i^l')) (70) 

SrTS = -^^{d.fW{d-{WTl')) (71) 
dp 

Kt - -^^{d.DMd^rT) (72) 

+ i^|^('I^T*') ■ ('^T^) + i^|l^(9''(#T^)) . (9,(#T*')) 

1 dGv .u . 1 9-Dv . . .a^ 
+ 2^^ > + 2^(^^>)(^ J ) 

+ 2^^^ ■ + 2^^^-^^-^ ■ 
~^{d.J.)u-{d^f) . (73) 

The nucleon isoscalar and isovector currents read 

^•M = ^^M^ ^ (74) 

= #r7''a' , (75) 

respectively. We write pu'^ — j'' and the four- velocity u'^ is defined as (1 — 
v^)~^^^{l,\r) where v is the three-velocity vector (v = in the rest-frame of 
the nuclear system). In addition to the isoscalar- vector S^, isoscalar-scalar Us, 
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isovector-vector St^ and isovector-scalar Sts self-energies, the density dependence 
of the vertex functions produces the rearrangement contributions SrS, S^, SrTS 
and S'^j,. The rearrangement terms result from the variation of the vertex function- 
als with respect to the baryon fields in the density operator p (which coincides with 
the baryon density in the nuclear matter rest- frame). For a model with density 
dependent couplings, the inclusion of the rearrangement self-energies is essential 
for energy-momentum conservation d^T'^'^ — 0, and thermodynamical consistency 
P^Wp (p) ~ i Si=i (i-®- ^'^^ pressure equation derived from the thermody- 
namic definition and from the energy-momentum tensor) 1471 1501 . 

When applied to nuclear matter or ground-state properties of finite nuclei, the 
point-coupling model is understood to be used in the mean-field approximation. The 
ground state of a nucleus with A nucleons is the product of the lowest occupied 
single- nucleon self-consistent stationary solutions of the Dirac equation l|68|l . The 
ground state energy is the sum of the single-nucleon energies plus a functional of 
the scalar density, 

A 

ps = ^ TpkTpk , (76) 
fc=i 

and of the baryon density of the nucleons, 

A 

p=Y.^U''' ^^^^ 

k = l 

where the sums run over occupied positive-energy single-nucleon states \k) with 
wave functions tpk- 

The density dependence of the coupling strengths which determine the self- 
energies 16617311 ■ will be constrained by in-medium QCD sum rules and chiral pion 
dynamics, to be specified in detail in the following sections 



4.2 Nuclear matter, part II 



In translationally invariant infinite nuclear matter all terms involving the derivative 
couplings 16H drop out and the single-nucleon Dirac equation reads 



where the self-energies are 

= Gvf , 



T-Sti,)-iM + Ss + T- EtsW 







X'y — GtvJt 1 
Es = Gs(>^*') , 
Sts = Gts{^t9) . 



(78) 



(79) 



(80) 

(81) 
(82) 
(83) 
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The total isoscalar vector self-energy includes the rearrangement contributions E!^. 
In nuclear matter at rest the spatial components of the four-currents vanish, and 
the densities are calculated by taking expectation values 

Ps = {•Pl'P^l'P) , (84) 
p = {<P\^')'"'P\$) , (85) 

Ps3 = {$\^T3^\$) , (86) 
P3 = {$\^T3'r°^\$) , (87) 

where \(I>) is the nuclear matter ground state. The energy density £ and the pressure 
P are derived from the energy-momentum tensor T'"' as 



£ = £kin + £^kin - ^GspI - \GtsP% + ^Gvp^ + \Gtv pl , (88) 



P = EpPp + E„p„ - f - £ki„ + ^Gvp^ + ^Grvpi + ^Gspl + ^GtsP% 
IdGs 2 IdGv 3 , IdGrv 2 , 1 dGrs 2 

The baryon density is related to the Fermi momentum fc/ in the usual way, 

d^k=^jg-, (90) 



^' (27r)3 J 37r2 

\k\<ky 

where the index i = {p, n} refers to protons and neutrons, respectively. The corre- 
sponding scalar densities are determined by the self-consistency relation 



i 2 r 3 m; m* 

P" = / d k- 



(27r)3 j l].2 + (M*)2 27r2 

|fc|<fc} 

with the proton and neutron quasi-particle energies 



(91) 



Ei= ^{k)Y-r{MlY , (92) 

and the effective nucleon masses 

M; = M + Gs Ps + Gts Ps3 , (93) 
M* = M + Gs Ps - Gts Pss ■ (94) 

The kinetic contributions to the energies of the protons and neutrons in nuclear 
matter are calculated from 

4in = ^ / d'kVk^T{Mf)^=\mp, + M*pl] i = {p,n}. (95) 
|fc|<fc» 

Note that rearrangement contributions appear explicitly in the expression for the 
pressure, whereas in the energy density, such contributions are implicit through the 
specific density dependence of Gv{p)- 
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4.3 QCD constraints 

We proceed now with a central theme of this section: establishing connections be- 
tween the density-dependent point couplings in the Lagrangian II58|I and constraints 
from QCD. Two key features of low-energy, non perturbative QCD are at the ori- 
gin of this discussion: the presence of a non-trivial vacuum characterized by strong 
condensates and the important role of pionic fluctuations governing the low-energy, 
long wavelength dynamics according to the rules imposed by spontaneously broken 
chiral symmetry. Our basic conjecture is, consequently, that the nucleon isoscalar 
self-energies 17911 and 1821 arise primarily through changes in the quark conden- 
sate and in the quark density at flnite baryon density, together with chiral (pionic) 
fluctuations induced by one- and two-pion exchange interactions. 

A. In-medium QCD sum rules 

QCD sum rules combine dispersion relation methods with operator product ex- 
pansion techniques in the analysis of correlations functions. When applied to the 
nucleon, the starting point is a suitably chosen 3-quark operator Jn{x) which car- 
ries nucleon quantum numbers. The correlation function 



has a pole structure which determines the nucleon mass in the presence of the 
(non-perturbative) QCD vacuum. That same correlation function can be written 
in terms of the operator product expansion. 



with vacuum expectation values of local operators On constructed from quark and 
gluon fields and organised according to their (mass) dimensions. The operators of 
lowest non-trivial dimension, d = 4, are qq and G^i^C^". The next higher ones 
include four-quark operators and combinations of quark and gluon fields. 

Now, for sufficiently large spacelike = —q^ > 0, the coefficients C„{q'^), 
called Wilson coefficients, decrease as inverse powers of with each of the higher 
dimensional operators to which they are attached. One can then apply QCD per- 
turbation theory to calculate the C„. The next step is to match 11 {—Q^) with its 
dispersion relation representation and apply a differential operation, the so-called 
Borel transform, to improve convergence properties. When this matching is per- 
formed, the leading order result for the nucleon mass in vacuum is loffe's formula 



Applications of QCD sum at finite baryon density p proceed in an analogous 
way, except that now operators such as q^g also have non- vanishing ground state 
expectation values: {q^q)p = 3p/2. In addition, the vacuum condensates vary with 
increasing density. Thus in-medium QCD sum rules relate the changes of the scalar 
quark condensate and the quark density at finite baryon density, with the isoscalar 
scalar and vector self-energies of a nucleon in the nuclear medium. In leading order 
which should be valid at densities below and around saturated nuclear matter, the 




(96) 




(97) 



n 
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condensate part, of the scalar self-energy is expressed in terms of the density 

dependent chiral condensate as follows |5||7||S|: 

4"..-2jK«).-(«).i^-f m 

with ps = We have used Eq.lBOJl together with the GOR relation GIJ- The 

chiral vacuum condensate {qq)o as a measure of spontaneous chiral symmetry break- 
ing in QCD is given in Eg. (1121 . The difference between the vacuum condensate {qqjo 
and the one at finite density involves the nucleon sigma term, crjv ~ {N\mqqq\N} , 
to this order. The Borel mass scale Ab ~ IGeV roughly separates perturbative and 
non-perturbative domains in the QCD sum rule analysis. 

To the same order in the condensates with lowest dimension, the resulting time 
component of the isoscalar vector self-energy is 

It reflects the repulsive density-density correlations associated with the time com- 
ponent of the quark baryon current, q'y'^q. Note that, as pointed out in Ref.|n|, the 
ratio 

(100) 



Zjg _ (TJV ps 



^(O'l 4(m„ + ma) p 

is approximately equal to —1 for typical values of the nucleon sigma term crjv and 
the current quark masses mu,d, and around nuclear matter saturation density where 
Ps — p (as an example, take crjv — 50 MeV and rriu + rad — 12 MeV at a scale of 1 
GeV). 

Identifying the free nucleon mass at p = according to loffe's formula II13II . 
M = ~^{qq)o, one finds 



Ef\p) ^M*{p)-M^ -^^Ps (101) 

and 



r^M- ^ p. (102) 



^(0)^^^ _ 4(m„ -I- md)M 

with the quark masses rriu^d to be taken at a renormalization scale p ~ Ab — 
inU =^ 1 GeV. 

Given these self-energies arising from the condensate background, the corre- 
sponding equivalent point-coupling strengths Gg^^ are simply identified through 
the relations 

4»)=Gfp. (103) 

and 

4"'=GrP (104) 

in mean field approximation. To this order, the couplings G^g y are density inde- 
pendent. Eg. llOlll implies (identifying M with the free nucleon mass): 



4°) ^ -350 MeV ^ . (105) 

50 MeV po 
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Evidently, the leading-order in-medium change of the chiral condensate is a source 
of a strong, attractive scalar field which acts on the nucleon in such a way as to 
reduce its mass in nuclear matter by more than 1/3 of its vacuum value as already 
anticipated in Ea. ll42ll . Using 11U3II one estimates: 



A correspondingly smaller value of Gg would result if only a fraction of the nucleon 
mass M is associated with the chiral condensate, leaving room for non-leading 
contributions from higher dimensional condensates. The QCD sum rules suggest 
that the ratio of the condensate scalar and vector self-energies is close to —1, so 
one expects roughly Gy ^ ~ — Gg^ 

One should note of course that the QCD sum rule constraints implied by 
Eas. lllUllllU21 and by the ratio IIUUII are not very accurate at a quantitative level. 
The leading-order loffe formula on which Eg. llOU relies has corrections from con- 
densates of higher dimension which are not well under control. The estimated error 
in the ratio S^^^ /Ey^ ~ — 1 is about 20%, given the uncertainties in the values of 
o"jv and rriu + md. Nevertheless, the constraints implied by Ea. mOOII give important 
hints for further orientation. 

B. Pionic fluctuations: self-energies from in- medium chiral 
perturbation theory 

The second essential ingredient is the contribution to the nucleon self-energies from 
chiral fluctuations related to pion-exchange processes (primarily two-pion exchange 
interactions, with small corrections from one-pion-exchange Fock terms). The cor- 
responding density dependent point-coupling strengths, denoted by G[^y{p) etc., 
are calculated using in-medium chiral perturbation theory.. 

In Section 3 we have shown how the nuclear matter equation of state can be 
calculated using in-medium chiral perturbation theory up to three loop order in the 
energy density, expanded in powers of the Fermi momentum kf (modulo functions 
of kf/m-n). The empirical saturation point, the nuclear matter incompressibility, 
and the asymmetry energy at saturation can be well reproduced at order ©(fc^) 
in the chiral expansion with just one single high-momentum cutoff scale which 
parametrizes short-distance physics. In this approach the nuclear matter saturation 
mechanism is entirely determined by the interplay between attraction from two-pion 
exchange processes and the stabilizing effects of the Pauli principle. 

The density dependence of the strength parameters G''^' (p) originating from pi- 
onic (chiral) fluctuations is calculated by equating the isoscalar-scalar, the isoscalar- 
vector, the isovector-scalar, and the isovector-vector self-energies II79I83II in the 
single-nucleon Dirac equation II78II with those calculated using in-medium chiral 
perturbation theory: 



4°) ~ -11 fm^ 



50 MeV 



at ps ~ po — 0.16 fm 



-3 



(106) 



(jg Ps 




(107) 
(108) 
(109) 
(110) 




Exploring the nucleus in the context of low-energy QCD 



31 



with 

(^) 1 dGg 2 1 dG^ 2 1 dGrs 2 1 i9GtV 2 

The equations of state of symmetric and asymmetric nuclear matter calculated in 
CHPT give, via the Hugenholtz - van Hove theorem, the sums U(T){kf , kf) — 
'^{T)s^(^/i P) + ^cryv^ikfjp) of the scalar and vector nucleon self-energies in 
the isoscalar and isovector channels at the Fermi surface p — kf up to two- 
loop order, generated by chiral one- and two-pion exchange IlOf. The differences 
^(T)s'^ikf,p) — ^{Tyv^ikf, p) are calculated from the same pion-exchange diagrams 
via charge conjugation (formally evaluating the anti-nucleon single particle poten- 
tial in nuclear matter). Following a procedure similar to the determination of the 
nucleon-meson vertices of relativistic mean-field models from Dirac-Brueckner cal- 
culations |49|, we neglect the momentum dependence of S'^r^^'^{p, p) and take their 
values at the Fermi surface p = kf. A polynomial fit up to order k^ is performed, 
and all four CHPT self-energies have the same functional form. 



C30 + C31A -f C32A^ -I- C3L In 



M2 



C50 + C5L In 



47r/^A 



(112) 



where the dimensionless parameter A is related to the momentum cut-off scale by 
A = 27r/^A, and /,r = 92.4 MeV denotes the pion decay constant. The cut-off 
parameter A ~ 0.65 GeV that has been adjusted to the nuclear matter saturation 
in Ref.|3], corresponds to A = 1.113. In order to determine the density dependence 
of the couplings from Eqs. (I107I11UII . the CHPT self-energies are re-expressed in 
terms of the baryon density p = 2fc^/37r^: 

r?«P^(fc/,p) = (c,l+C,2/5* +Cs-ipi)p, (113) 

r{?"P^(fc/,p) = (C„i +C„2P* +C„3P*)p, (114) 

S?S^'^{kf,p) = {Ctsl+Cts2pi +CtsSpi)P3 , (115) 

SxW^ikfyP) = {Ctvl + Ctv2P^ +Cti,3p3)p3 . (116) 

The resulting expressions for the density dependent couplings of the pionic fluctu- 
ation terms are* 

G'-g^ = Csl + Cs2P^ + CsSP^ , (117) 



i2 — Cs2 — I 
C„3 = (3c„3 - Cs3 - (5^(cts3 + Ct„3))/4 



G^' = Cvl + Cv2p'> + Cv3p^ < C„2 = (6c„2 " Cs2 — 5^(cts2 + Ctv2))/7 ,(118) 



= + '=*''2P3 + Cts3P^ , (119) 

G!j!1 = Ctvl + Ctv2P^ + Ct„3P* , (120) 

where 5 — {pp — Pn)/{pp + Pn)- For A = 1.113, the coefficients of the expansion of 
the CHPT self-energies in powers of the baryon density Eas. 1113111611 are given in 
Table [U 

* small differences between ps and p at nuclear matter densities are neglected here 
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Table 1. The coefficients of the expansion of the in-medium CHPT self energies 
111311 - 111611 in powers of the baryon density, for the value of the cut-off parameter 
A = 646.3 MeV. In this case we assumed N = Z, neglecting isovector contributions. 



Csl 


-2.805 fm^ 


Ctsl 


-0.373 fm^ 


Cs2 


2.738 fm^ 


Cts2 


0.251 fm^ 


Cs3 


1.346 fm* 


Cts3 


0.387 fm'' 


Cvl 


-2.718 fm^ 


Ctul 


-0.562 fm^ 


Cv2 


2.841 fm^ 


Ctv2 


0.207 fm^ 


Cv3 


1.325 fm^ 


Ctv3 


0.398 fm'' 



Combining now the effects from leading order QCD condensates and pionic 
fluctuations, the strength parameters of the isoscalar four-fermion interaction terms 
in the Lagrangian 15811 are: 

Gs.v(p) = 4%+4:i-(p) • (121) 

For the isovector channels we assume that only pionic (chiral) fluctuations con- 
tribute. 

C. Corrections of higher order 

Up to this point the nucleon self-energies are evaluated at chiral two-loop level 
including terms of order k^. At order (oc p^), several additional effects appear. 
First, condensates of higher dimension enter, such as {qFqqrq), {q^q)^. Their de- 
tailed in-medium dependence is difficult to estimate. Secondly, four-loop CHPT 
contributions to the energy density (three-loop in the self-energies) introduce gen- 
uine 3-body interactions. All those effects combine to produce the C'(fc/) correction 
in the energy per particle. Our conjecture is that these terms should be subleading 
in the fc/-expansion. In estimating those contributions, we will take a pragmatic 
point of view. Given the density-independent leading condensate terms G*-"^ and 
the CHPT couplings G^^'^p) evaluated to 0{k)) ix p3 ^ we generalize 

G(p) = G'"' -I- G<"' (p) + <5G('' (p) , (122) 

adding terms SG^^^ = g^^^ p and let the constant p'^' be determined by a least 
squares fit to properties of nuclear matter and finite nuclei. At first sight, the cor- 
rections dG^g\ 5Gy\ stc... would then introduce several additional parameters, 
undermining our startegy to keep the number of freely adjustable fine-tuning con- 
stants as small as possible. For the scalar self-energy a primary uncertainty arises 
from the four-quark condensate. It is common practice to approximate (qqqq) as- 
suming factorization into the form {qq)'^, which introduces potentially large and 
uncontrolled errors as discussed in detail in Ref. jgj . While the in-medium values of 
condensates such as {qFqqFq) are thus not well determined, it nevertheless appears 
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that only a weak or moderate density dependence of these condensates is consistent 
with known nuclear phenomenology. We do not have to be much concerned about 
this issue here because our explicit treatment of scalar tttt fluctuations (up to order 
in the energy density) removes at least part of such uncertainties. What we can 
conclude from this calculation is that factorization of the four-quark condensates, 
assuming ground state dominanance, is not justified since the QCD ground state 
strongly couples to the scalar two-pion continuum via the operator qq. 

So as a first option, we assume that the only non-vanishing ^G'^-* is SGy^ = 
gy ^ p, and determine the constant gy \ Note that this SG'^'^^ produces a self-energy 
contribution SE^^^ oc which acts like a three-body force. 

4.4 Nuclear matter, part III 

We return now to the calculation of the nuclear matter equation of state (EOS), 
using the relativistic point-coupling model with QCD-constrained input as outlined 
in the previous sections. These calculations are performed in three steps. 

i) In the first step we follow the implicit assumption made in Ref.|2] and in 
Section 3, namely that X"^"' = — X"^' in nuclear matter, and ignore the contribu- 
tion of the condensate background self-energies. The point-coupling model, with 
the density dependence of the strength parameters determined by pionic (chiral) 
fluctuations, nicely reproduces the nuclear matter EOS calculated previously us- 
ing in-medium CHPT: the binding energy per particle, the saturation density, the 
compression modulus, and the asymmetry energy. Small differences arise mainly be- 
cause in the mapping of the CHPT nucleon self-energies on the self-energies of the 
point coupling model in Eqs. llUTIlim , the momentum dependence of the former 
has been frozen to the values at the Fermi momentum. This is a well known prob- 
lem that arises also in the determination of the meson-nucleon in-medium vertices 
from Dirac-Brueckner calculations of nucleon self-energies in nuclear matter [49 1 . 

At very low-densities, the CHPT equation of state cannot reproduce realistic 
calculations of E/A, for obvious reasons. The low density limit of the energy per 
particle is not given by a free Fermi gas of nucleons. Proton-neutron pairs coalesce to 
form deuteron bound-states, an entirely non-perturbative phenomenon that cannot 
be handled at any order of a perturbative approach. The low-density limit of E/A 
approaches a constant as fc/ ^ 0, namely one half of the deuteron binding energy. 
The effect of this shift is marginal, however, for the bulk of nuclear matter around 
saturation density. 

ii) The next step includes the contributions of the condensate background self- 
energies, S^gy, in the isoscalar-scalar and isoscalar- vector channels. In Ref. )12l we 
have shown that, even though to first approximation the condensate potentials 
do not play a role in the saturation mechanism, they are indeed essential for the 
description of ground-state properties of finite nuclei. 

At this stage we first treat both Gg"' and Gy' as independent parameters, 
irrespective of the QCD sum rule constraint IIUUII . The three parameters: Gg"', 
Gy^ and A (the CHPT cutoff scale), are then adjusted to reproduce "empirical" 
nuclear matter properties: E/A = -16 MeV (5%), po = 0.153 fm"^ (10%), the 
compression modulus Ko = 250 MeV (10%), and the asymmetry energy Ao — 33 
MeV (10%). The values in parentheses refer to the error bars used in the fitting 
procedure. 



34 Dario Vretenar and Wolfram Weise 



Table 2. Nuclear matter saturation properties calculated in the relativistic point- 
coupling model constrained by in-medium QCD sum rules and chiral perturbation 
theory. In addition to the chiral one- and two-pion exchange contribution to the 
density dependence of the coupling parameters, the nuclear matter EOS shown in 
the first row incorporates the isoscalar condensate background nucleon self-energies 
linear in the corresponding densities. The EOS displayed in the second row is cal- 
culated by including also the non-linear "3-body" contribution (proportional to p^) 
in the isoscalar vector condensate self-energy. 



E/A[MeV] po[fm-^] A'o[MeV] M*/M ylo[MeV] 



LINEAR -15.97 0.148 283 0.753 45.3 
NON-LINEAR -15.76 0.151 332 0.620 30.2 



The corresponding least-squares /it yields G^'' = —7 fm^, G^' = 7 fm^, and A = 
685 MeV. The result is shown in Fig. |1] (EOS of isospin symmetric nuclear matter 
in comparison with the one calculated exclusively from chiral pion-exchange) and in 
the first row of Table|5| The EOS so obtained is not yet satisfactory. The asymmetry 
energy Aq at saturation is too high, and the values of G^"' and G^' are smaller 
than the leading order QCD sum rule estimate 1)106^ . This results in the relatively 
large Dirac effective mass M* /M — 0.75. A high Dirac mass in turn indicates 
that the effective single-nucleon spin-orbit potential is too weak to reproduce the 
empirical energy spacings between spin-orbit partner states in finite nuclei. Stronger 
background scalar and vector fields are required in order to drive the large spin-orbit 
splittings in finite nuclei. 

Notice however that, even though Gf^ and G^y^ were varied independently, the 
minimization procedure chooses to balance the contributions from the correspond- 
ing large scalar and vector self-energies, such that their sum tends to vanish: there 
is already sufficient binding from pionic fluctuations alone, so that even the uncon- 
strained fit prefers S^g ^ = —Sy ^ for the condensate potentials. This remarkable 
feature prevails as we now proceed with fine-tuning improvements. 

iii) The third step incorporates in addition the higher-order correction SGy ^ = 
g'^^p, see Eg.lin^. For Gf^ = -12 fm^ G^' = 11 fm^ g'^^ = -3.9 fm^ and 
A = 600 MeV, determined in a least-squares fit to the "empirical" nuclear matter 
input, the EOS resulting is displayed in Fig. 01 and the nuclear matter properties 
at saturation are listed in the second row of Table |5| This nuclear matter EOS is 
actually quite satisfactory, with realistic values of the binding energy, saturation 
density, asymmetry energy at saturation, and a low effective Dirac mass. The only 
exception is a relatively large nuclear matter incompressibility. It could be cured 
by introducing a p"^ term in the expansion of the self-energies, but we prefer not to 
invoke additional fine-tuning parameters at this level of the discussion. 




Fig. 4. Binding energy per nucleon for symmetric nuclear matter as a function of 
the baryon density, calculated from chiral one- and two-pion exchange between nu- 
cleons (case 1, solid curve), by adding the isoscalar condensate background nucleon 
self-energies linear in the corresponding densities, with G^^' = — Gg""* — 7 fm^ (case 
2, dotted curve), and finally by including also the non-linear contribution g^'p^ to 
the isoscalar vector condensate self-energy, with Gy^ — 11 fm"^, Gg"' = —12 fm^ 
(case 3, dashed curve). 



The resulting couplings Gg y of the condensate background fields are remark- 
ably close to the prediction lilU6ll of the leading order in-medium QCD sum 
rules. The scalar and the vector "condensate" self-energies S^g-' = Gg^'ps and 
Sy^ = Gy^p follow the expectation E^^^ ~ —Ey'^ to within 5% at saturation 
density, even without this condition being pre-imposed. The large isoscalar conden- 
sate background self-energies in turn lead to a relatively low effective Dirac mass, 
crucial for the empirical spin-orbit splittings in finite nuclei. Finally, the cut-off 
A = 600 MeV differs by less than 10 % from the value (646 MeV) obtained when 
the nuclear matter EOS results exclusively from one- and two-pion exchange be- 
tween nucleons The &Gy^ — gy^ p term is relatively small: at saturation density, 
5Gy^ /Gy^ ~ 0.05. Splitting JG'^' between scalar and vector parts (e.g. by choosing 
5Gy^ ~ ^Gg^' ~ —2 fm~^) would induce only marginal changes. 
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4.5 Comparison with Dirac-Brueckner G-matrix theory 

At this point it is instructive to compare the density-dependence of our self-energies 
Ss,v{p) with that of the scalar and the vector self-energies resulting from Dirac- 
Brueckner G-matrix calculations which start from realistic NN- interactions VArjv- 
Brueckner calculations iterate Vjvjv to all orders in the ladder approximation, keep- 
ing track of Pauli exclusion priciple effects on intermediate A'^A'^ states. In addition, 
the Dirac-Brueckner approach takes into account the lower components of the rel- 
ativistic (Dirac) wave functions of the nucleons, so that the density dependent 
nucleon self-energies resulting from such calculations separate into Lorentz scalar 
and vector parts. 




Fig. 5. Comparison of the fc/-dependencies of isoscalar vector and scalar self- 
energies resulting from a Dirac-Brueckner calculation (solid lines: DB 52 ) with 
the self-energies generated from in-medium chiral perturbation theory (dashed lines: 
CHPT -I- point-couplings) up to 3-loop order in the energy density. (We thank Paolo 
Finelli for preparing this figure.) 



We refer here explicitly to a calculation, reported in Ref. which uses the 
Bonn A potential. With this potential, the nuclear matter saturation density comes 
out somewhat too high {psat — 0.185 fm~'^), higher than the one of our best fit. In 



Exploring the nucleus in the context of low-energy QCD 



37 



order to perform a meaningful comparison, simulating this larger saturation den- 
sity in our approach requires weakening the condensate mean fields X'g"^ by about 

25 %, but with no changes in the pionic terms S^^y. After this readjustment the 
difference in the fc/-dependences between our Ss and and those resulting from 
the Dirac-Bruckner calculation is less than 10 % over the entire range of densi- 
ties from 0.5po to 2.5po (see Fig. |^. This is a remarkable observation: it appears 
that in-medium chiral perturbation theory at two-loop order, with a cut-off scale 
A ~ 0.6 GeV that converts unresolved short-distance dynamics at momenta be- 
yond A effectively into contact terms, generates quantitatively similar in-medium 
nucleon self energies as a full Dirac-Brueckner calculation, when requiring that both 
approaches reproduce the same nuclear matter saturation point. The reasoning be- 
hind this observation can presumably be traced back to the separation of scales at 
work, a key element of chiral effective field theory on which our approach is based. 
One- and two-pion exchange in-medium dynamics at momentum scales compara- 
ble to kf must be treated explicitly. Our in- medium CHPT approach includes all 
terms (ladders and others) to three-loop order in the energy density, a prominent 
one being iterated one-pion exchange. The non-trivial kf dependence in the self- 
energies (beyond "trivial" order fc^) reflects the action of the Pauli principle in 
these processes. Such features are also present in leading orders of the Brueckner 
ladder summation which includes, for example, Pauli effects on iterated one-pion 
exchange. On the other hand, the iteration of short-distance interactions to all or- 
ders involves intermediate momenta much larger than kf and generates effectively 
the same self-energy pieces (proportional to density p) as a suitably chosen con- 
tact interaction: the infinite ladder summation of short-distance interactions can 
simply be subsumed in a contact term reflecting the high-momentum scale A. The 
reasoning here has a close correspondence to related studies in Ref. |53|. 



4.6 Intermediate summary and discussion 

While the quantitative details for the terms contributing to the scalar and vector 
self-energies can be reconstructed from the tables and figures that we have included 
in this section, it is instructive to observe systematic trends in these self-energies 
when expanded in powers of kf (or equivalently, in powers of p3). We focus on 
symmetric nuclear matter and write 

Es,v = + ^s,l + S^s,v , (123) 

combining the leading condensate terms S^s\^, the chiral (pionic) terms s]^y cal- 
culated to order fcp and possible corrections of higher order summarized in SSs.v- 
We introduce the saturation density po ~ 0.16 fm~^ as a convenient scale and write 
the condensate terms as 

4«)(p).-0.35Gevif|i (g), (124) 

r(?)(p).-,0.35Gevif|l (ii). (125) 

We note that with |Gg°''| chosen less than 5% larger than Gy', this just compen- 
sates for the difference between the baryon density p and the scalar density ps < p 
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such that X'g"-' 



-,(0) 



results at p = po- The "best-fit" values Ci"'' 



^12 fm^ 



and G 



(0) 



-l-llfm are thus perfectly consistent with QCD sum rule expecta- 



tions: a non-trivial result. The pionic fluctuations from one- and two-pion exchange 
processes have the following approximate pattern: 



-75 MeV 



1 + di 



+ d2 



where di varies between —0.61 (scalar) and —0.65 (vector) and d2 
cases. Finally, the higher order corrections are summarized as 



(126) 



-0.17 for both 



6Es + 5Ev 



-20 MeV 



(127) 



where the detailed decomposition into scalar and vector part turns out not to be 
relevant. 

The leading attraction at 0{kj) in E^^y depends on the cut-off scale A = 27r/,r A 
which separates "active" two-pion exchange dynamics at long and intermediate dis- 
tances from short-distance (high momentum) dynamics. Details at short distance 
scales (related to the intrinsic structure of the nucleons and their short-range in- 
teractions) are not resolved as long as the Fermi momentum fc/ is small compared 
to the chiral symmetry breaking scale, 47r/,r. The unresolved short-distance infor- 
mation can then be translated into an equivalent (density-independent) four-point 
vertex which generates self-energy terms linear in the density at the mean-field 
level. 



5 Finite Nuclei 

In this section the relativistic nuclear point-coupling model, constrained by in- 
medium QCD sum rules and chiral perturbation theory, is applied in calculations of 
ground state properties of finite nuclei. In the mean-field approximation the ground 
state of a nucleus with A nucleons is represented by the antisymmetrized product of 
the lowest occupied single- nucleon stationary solutions of the Dirac equation 16311 . 
The calculation is self-consistent in the sense that the nucleon self-energies are func- 
tions of nucleon densities and currents, calculated from the solutions of Eq. 16311 . 
The mean-field approach to nuclear structure represents an approximate implemen- 
tation of Kohn-Sham density functional theory (DFT) 54 55 56 , which is suc- 
cessfully employed in the treatment of the quantum many-body problem in atomic, 
molecular and condensed matter physics. At the basis of the DFT approach are en- 
ergy functionals of the ground-state density. In relativistic mean-field models these 
become functionals of the ground-state scalar density and of the baryon current. 
The scalar and vector self-energies play the role of local relativistic Kohn-Sham 
potentials ^[5]. The mean-field models approximate the exact energy functional, 
which includes all higher-order correlations, by powers and gradients of fields or 
densities, with the truncation determined by power counting [5]. 
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5.1 Ground state energy 

In this lecture wc only consider spherical even-even nuclei. Because of time- reversal 
invariance the spatial components of the four-currents vanish, and the nucleon self- 
energies (reduced to the time component of the vector terms S'^) read: 

i:° = GvP + DvAp-eA°^^^-^ , (128) 



Gtv P3 + DrvApa , (129) 

Ss = GsPs + DsAps , (130) 
StS = GtS Ps3 + DtsAps3 , (131) 

^.s = ^(Vp,).(Vp), (132) 

SrTS = ^^{Vps3)-iVp), (133) 

S°T = ^^i^P^) ■ (Vp) , (134) 

, 1 dGs 2 1 dDs . ,^ . 
''^=+2^^=-2^(^^^)-(^^^) 

+ 2^^ -2^(^^)-(^^) 

1 OGts 2 1 ODts 
+ 2^^=^-2^(^^^^)-(^^^^) 

1 OGtv 2 1 dDrv ^ ^ 
+ 2^''^-2^(^^^)-(^^^) 

+^(Vp)-(Vp). (135) 

In the mean-field approximation the isoscalar-scalar, isoscalar- vector, isovector- 
scalar and isovector-vector ground state densities are calculated from the Dirac 
wave functions ipk of the occupied positive-energy orbits as: 

A 

= ^ V'fcV'fc , (136) 
fc=i 

A 

P = ^4i^k, (137) 
fc=i 

A 

Ps3 = ^ '>PkT3lpk , (138) 
fc=l 

A 

P3 =ZlV'lT-3Vfe , (139) 



respectively. The expression for the ground state energy of a nucleus with A nucleons 
reads: 
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A „ 

£ = ^ tfc - i / cPx{Gspi + Ds PsAps + Gvp" + DvpAp 

+GtvpI + DrvpzApz + Gts plz + Dts Ps3 Ap^g + fVc 
dGs 2 , dGv 3 , dGrv 2 , dGrs 2 

■ (VP.)P+ ^(Vp) • (Vp)p 

+2^(Vp.) ■ (Vp)p. + 2^^(Vp3) ■ (Vp)p3 

+ 2^|^(Vp.3)-(Vp)p.3| , (140) 

where ek denotes the single-nucleon energies. In addition, for open shell nuclei pair- 
ing correlations are included in the simple BCS approximation jFTj . After the solu- 
tion of the self-consistent Dirac equation, the microscopic estimate for the center- 
of-mass correction 

is subtracted from the total binding energy. Here Pcm is the total momentum of 
the nucleus with A nucleons. 



5.2 Surface (derivative) terms 

The calculated properties of finite nuclei depend of course on an accurate tuning 
of coupling parameters. We keep the variational freedom in those parameters at 
minimum, constraining the density-dependent couplings Gs, Gv, Gts and Gtv 
by in-medium QCD sum rules and explicit chiral perturbation theory calculations 
of one- and two-pion exchange diagrams, as described in the previous sections. For 
finite nuclei, fine-tuning of their detailed surface structure will be required, so we 
must also determine the coupling parameters of the derivative terms: Ds, Dy, Dts 
and Dtv ■ Dimensional considerations suggest the following ansatz 

Dip) = ^ , (142) 

where M is a. characteristic mass scale for a given spin-isospin channel. There is 
no deeper reason, however, for the derivative terms to have the same density de- 
pendence as the coupling parameters of the four-fermion interactions. The density 
dependence of Di{p) could, in principle, be derived by CHPT calculations for in- 
homogeneous nuclear matter. A simpler option, followed here, is to treat the D's 
as density-independent adjustable constants. In this case the remaining rearrange- 
ment contribution to the vector-self energy becomes, of course, much simpler. As 
it has been emphasized by Serot and Furnstahl the empirical data set of bulk 
and single-particle properties of finite nuclei can only constrain six or seven param- 
eters in the general expansion of the effective Lagrangian in powers of the fields 
and their derivatives. In particular, only one parameter of the derivative terms can 
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Fig. 6. Neutron and proton single-particle levels in O and Ca calculated in the 
relativistic point-coupling model. The density dependent coupling strengths include 
only the contribution from chiral one- and two-pion exchange between nucleons. 
(From Ref. [T^.l 



be determined by the binding energies and radii of spherical nuclei. In the present 
analysis we therefore set Dv, Dts and Dtv equal to zero, and adjust the single 
remaining surface parameter Ds of the isoscalar-scalar derivative term to prop- 
erties of light and medium-heavy N ~ Z nuclei. This approximation, which was 
first used by Serot and Walecka in Ref. |58| . implies that the isoscalar- vector, the 
isovector-scalar and the isovector-vector interactions are considered to be purely 
contact interactions (no gradient terms). 



5.3 Single-particle energies 

It is instructive to consider separately the contributions from chiral pion dynamics 
and condensate background self-energies to properties of finite nuclei. In a first 
step we have calculated the ground states of and '^^C& using the coupling 
parameters determined by the nuclear matter EOS of Ref. Gs{p) = G^g\p), 
Gv{p) = G\;\p), Gts{p) = G^''i(p), Gtv{p) = G^"^-(p), A = 646.3 MeV, while 
the couplings Gg^y to the condensate background fields are set to zero. In this case 
the nuclear dynamics is completely determined by chiral (pionic) fluctuations. The 
resulting total binding energies are already within 5—8 % of the experimental values, 
but the radii of the two nuclei are too small (by about 0.2 fm). This is because the 
spin-orbit partners (IP3/2, lPi/2) and (1^5/2,1^3/2) are practically degenerate. In 
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Fig.|S]we display the calculated neutron and proton single-particle levels of and 
''"Ca. The energies of the degenerate doublets are close to the empirical positions of 
the centroids of the spin-orbit partner levels, and even the calculated energies of the 
s-states are realistic. This is an interesting result. Chiral pion dynamics provides 
the saturation mechanism and binding of nuclear matter, but not the strong spin- 
orbit force. The inclusion of the isoscalar-scalar derivative term has some effect on 
the calculated radii but it cannot remove the degeneracy of the spin-orbit doublets. 
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Fig. 7. The neutron and proton single-particle levels in O calculated in the 
relativistic point-coupling model, are shown in comparison with experimental levels. 
The calculation is performed by including both the contributions of chiral pion- 
nucleon exchange and of the isoscalar condensate self-energies. (From Ref. 



The spin-orbit potential plays a central role in nuclear structure. It is at the 
basis of the nuclear shell model, and its inclusion is essential in order to reproduce 
the experimentally established magic numbers. In non-relativistic models based 
on the mean-field approximation, the spin-orbit potential is included in a purely 
phenomenological way, introducing the strength of the spin-orbit interaction as an 
additional parameter. Its value is usually adjusted to the experimental spin-orbit 
splittings in spherical nuclei, for example ^^O. On the other hand, in the relativistic 
description of the nuclear many-body problem, the spin-orbit interaction arises 
naturally from the scalar-vector Lorentz structure of the effective Lagrangian. In 
the first order approximation, and assuming spherical symmetry, the spin-orbit 
term of the effective single-nucleon potential can be written as 
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(143) 



where the large spin-orbit potential Vis arises from the difference of the vector and 
scalar potentials, V - S ^ 0.7 GeV 01121. Explicitly, 



Vis 



M 



M, 



eff 



(V^S) 



(144) 



where M^ff is an effective mass specified as |5] 



M,ff = M~-{V~S). 



(145) 



The isoscalar nucleon self-energies generated by pion exchange are not sufficently 
large to produce the empirical effective spin-orbit potential. The degeneracy of spin- 
orbit doublets is removed by including the isoscalar background self-energies S'^'y 
that arise through changes in the quark condensate and the quark density at finite 
baryon density. 
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Fig. 8. Same as in Fig.EI but for ■*''Ca.(From Ref. P).! 



The effect of including S^ y is demonstrated in Figs. |7| and |S| In addition to 
the four parameters (G^°' = -12 fm^ G*^' = 11 fm^ g'^^ = -3.9 fm^ and A = 600 
MeV) determined by the nuclear matter equation of state, we have adjusted the 
isoscalar-scalar derivative term: Ds — —0.713 fm*. This value of Ds is very close to 
the ones used in the effective interactions of the standard relativistic point-coupling 
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Fig. 9. Neutron and proton single-particle spectra of ^^Ni, calculated with the 
standard relativistic mean-field model NL3 effective interaction |61|. and with the 
relativistic point-coupling model constrained by in-medium QCD sum rules and 
chiral perturbation theory (FKVW [SI). 



model of Ref. |41|. It is also consistent with the "natural" order of magnitude 
expected from Ds ~ Gs/A^. 

The calculated neutron and proton single-particle energies of and ^''Ca are 
in excellent agreement with the empirical single-nucleon levels in the vicinity of 
the Fermi surface. In particular, with the inclusion of the isoscalar condensate self- 
energies, the model reproduces the empirical energy differences between spin-orbit 
partner states. This is an important result. It supports our primary conjecture: 
while nuclear binding and saturation are almost completely generated by chiral 
(two-pion exchange) fluctuations in our approach, strong scalar and vector fields 
of equal magnitude and opposite sign, induced by changes of the QCD vacuum in 
the presence of baryonic matter, generate the large effective spin-orbit potential in 
finite nuclei. Not surprisingly, the '*''Ca spectrum is reminiscent of an underlying 
pseudo-spin symmetry |60l . 

In Fig. |5]we compare the single-nucleon spectra of ^^Ni, calculated in our ap- 
proach, with the results of a relativistic mean-field calculation using NL3 (61,|, prob- 
ably the best phenomenological non-linear meson-exchange effective interaction. 
The agreement is convincing. All these results demonstrate that in the present ap- 
proach, based on QCD sum rules and in-medium chiral perturbation theory, and 
with a small number of model parameters determined directly by these constraints, 
it is possible to describe symmetric and asymmetric nuclear matter, as well as prop- 
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erties of finite nuclei, at the same quantitative level as the best phenomenological 
relativistic mean-field models. 

5.4 Systematics of binding energies and charge radii 

In Table El we summarize our calculated binding energies per nucleon and charge 
radii of light and medium-heavy nuclei, in comparison with experimental values. 
The parameters have been kept unchanged from those used in the ^^O and ^''Ca 
calculations. The resulting agreement between the calculated and empirical binding 
energies and charge radii is indeed very good. 

Table 3. Binding energies per nucleon and charge radii of light and medium-heavy 
nuclei, calculated in the relativistic point-coupling model constrained by in-medium 
QCD sum rules and chiral perturbation theory are compared with experimental 
values. 





£/A''^P (MeV) E/A (MeV) 


r^^P (fm-^) 


Tc (fm-^) 




7.976 


8.027 


2.730 


2.735 


^°Ca 


8.551 


8.508 


3.485 


3.470 


"^Ca 


8.617 


8.537 


3.513 


3.473 


*«Ca 


8.666 


8.964 


3.484 


3.486 




8.260 


8.182 




3.551 




8.756 


8.779 




3.571 


s^Cr 


8.776 


8.635 


3.647 


3.641 


58 Ni 


8.732 


8.493 


3.783 


3.778 




8.777 


8.775 


3.868 


3.879 


««Sr 


8.733 


8.855 


4.206 


4.234 


90 Zr 


8.710 


8.746 


4.272 


4.284 



5.5 Density distributions 

Calculated nucleon densities of ^^O and ^^Ca are plotted in Figs. [TOl and [TT] The 
results obtained in our QCD-constrained relativistic point-coupling model |13| are 
compared with the nucleon densities calculated in the standard relativistic mean- 
field meson- exchange model with the NL3 effective interaction |61|. Also shown are 
the corresponding charge densities in comparison with the empirical charge density 
profiles |62| . When compared not only with NL3, but also with other standard 
relativistic mean-field effective interactions, the density profiles calculated in our 
point-coupling model display less pronounced shell effects, in better agreement with 
empirical densities. 
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Fig. 10. Nucleon density of O as a function of tlie radial coordinate. The result 
obtained in the QCD-constrained relativistic point-coupling model (FKVW 
is compared with the nucleon density calculated in the standard relativistic mean- 
field model with the NL3 effective interaction |61) . In the insert the corresponding 
charge densities are compared with the empirical charge density profile |62| . 

6 Concluding remarks and outlook 

Our aim in these lectures notes has been to demonstrate the key role of a guid- 
ing principle, emerging from the non-perturbative dynamics of QCD: spontaneous 
chiral symmetry breaking. The effective field theory based on this principle ties 
together a wide variety of strong-interaction physics involving the lightest quarks: 
from the low-energy interactions of pions and nucleons via the condensate structure 
of the QCD vacuum and its global thermodynamical implications, to fundamental 
aspects of nuclear binding and saturation. Much work lies ahead; the exploration of 
the nuclear many-body problem in view of its low-energy QCD connections promises 
new insights. At the least, the QCD constraints can significantly reduce the freedom 
in the choice of parameters. The wide variety of results reported in previous sections 
have been obtained with only five parameters, two of which (the couplings C'^g y) 
turn out to be so remarkably close to leading order QCD sum rule estimates that 
one could have guessed their values right from the beginning. The remaining three 
parameters (the high-momentum scale refiecting unresolved short-distance physics. 
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Fig. 11. Same as in Fig.HSl but for ''°Ca. 

a "three-body" contribution to the nucleon self-energy and a surface (derivative) 
term) behave "naturally" according to the power-counting doctrine of effective field 
theory. 

Of course, questions about systematic convergence of the in-medium chiral loop 
expansion still remain and need to be further explored. In particular, the important 
role of the zi(1230), the prominent spin-isospin excitation of the nucleon, is so far 
just hidden in the parametrization of "short-distance" phenomena. However, the 
scale associated with the mass difference Ma — Mn ~ 0.3 GeV is in fact comparable 
to the nuclear Fermi momentum, suggesting a more explicit treatment of the A in 
the chiral effective Lagrangian. 

Calculations for heavier N ^ Z nuclei are in progress. Systems such as '^"^Pb 
and beyond are a test case for the detailed isospin structure of the underlying 
chiral pion dynamics. In the present version of the model, the isovector parts of the 
interaction are exclusively determined by chiral one- and two-pion exchange. Global 
properties of asymmetric nuclear matter and neutron matter are reproduced quite 
well at least at moderate densities, but pion exchange processes do not account 
for the short-range part of the isovector effective interaction: the resulting density 
dependent couplings Gts and Gtv are presumably still too weak. There is need for 
stronger self-energies in the isovector channel. In meson- exchange language, such 
additional short-distance interactions could arise from the p and ao resonances 
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(where we must note, of course, that a large non-resonant low-mass part of the p 
meson channel is already accounted for by explicit isovector two-pion exchange). 
This is an important topic, to be further addressed in future studies which will 
focus on a quantitative description of heavy nuclei and extrapolations into regions 
of extreme isospin. 
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